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Abstract Flattened, rotating disks of cool dust and gas extending for tens to hundreds of AU 
are found around almost all low mass stars shortly after their birth. These disks generally persist 
for several Myr, during which time some material accretes onto the star, some is lost through 
outflows and photoevaporation, and some condenses into centimeter- and larger-sized bodies 
or planetesimals. Through observations mainly at infrared through millimeter wavelengths, we 
can determine how common disks are at different ages, measure basic properties including mass, 
size, structure, and composition, and follow their varied evolutionary pathways. In this way, we 
see the first steps toward exoplanet formation and learn about the origins of the Solar System. 
This review addresses observations of the outer parts, beyond 1 AU, of protoplanetary disks 
with a focus on recent infrared and (sub-)millimeter results and an eye to the promise of new 
facilities in the immediate future. 
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1 INTRODUCTION 

Circumstellar disks are an inevitable consequence of angular momentum conser- 
vation during the formation of a star through gravitational collapse. Initially 
disks rapidly funnel material onto the star but, as the surrounding molecular 
core is used up or otherwise disperses, the accretion rate decreases and a small 
amount of material persists. That these disks can be considered protoplanetary 
is apparent not only through the geometry of the Solar System but also the high 
detection rate of exoplanets. 

Because disks exhibit a range of temperatures - hot near the star, cooler father 
away - they radiate strongly at a range of wavelengths from microns to millime- 
ters. They can therefore be observed with infrared and radio telescopes, and the 
mapping of wavelength to radius allows detailed models of their structure to be 
determined purely from unresolved photometry. Furthermore, their longevity, 
relative to the natal core, allows their properties to be studied in relation to the 
optically visible protostar. 

Internal friction, or viscosity, within the disk drives continued accretion onto 
the star. To preserve angular momentum, some material is lost through outflows 
and the disk may gradually spread out with time. Its structure may also be 
strongly affected by photoevaporation, both from the central star and external 
stars, and by the agglomeration of dust grains well beyond the typical sizes found 
in the interstellar medium including, ultimately, into planetesimals large enough 
to gravitationally perturb the disk. The various evolutionary pathways lead to 
inner holes and gaps that reveal themselves through a relative decrement in flux 
over a narrow range of wavelengths and which may also be imaged directly at 
sufficiently high resolution. 

The Infrared Astronomical Satellite (IRAS) opened up the infr a,red sky and 



allowed the first statistical studies of disk occurrence to be made ( Strom et al 



19891 ). Shortly thereafter the first sensitive detect ors at millimeter wavelengths 



show ed that many disks contained large dust grains (IWeintraub. Sandell Duncan 



19891) with enough m aterial to form planetary systems on the scale of our own 



Beckwith etld] Il99nl ). Interferometry at t hese long wavelengths pro vided the 



ability to resolve the rotation in the disks ( Sargent Sz Beckwith 198?! ). but un- 



equivocal evidence for their fiattened morphology actually came in the optical. 
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with the Hubble Space Telescope (Hubble), t hrough exquisite im ages of disk 
shadows against a bright nebular background ( O'dell &: Wen 19941 ). The pace 



of discoveries has accelerated in the past decade due to increases in sensitiv- 
ity, resolution, and wavelength coverage. The Infrared Space Observatory {ISO) 
and, in particular, the Spitzer Space Telescope (Spitzer) have greatly expanded 
the known disk inventory in terms of central stellar mass, age, environment, 
and evolutionary state. Interferometry has expanded to longer baselines and 
shorter wavelengths, including into the sub-millimeter regime with the Submil- 
limeter Array (SMA), providing the ability to map fainter structures in greater 
detail. The potential to address fundamental questions in protoplanetary disk 
studies provided significant motivation for the development of major new facil- 
ities including the Herschel Space Observatory {Herschel) and Atacama Large 
Millimeter /submillimeter Array [ALMA). 

The rather short history of protoplane tary disk research may be followed in the 
regular Protostars and Planets series ( Black & : Matthew ^ ^1985. Gehrels 1978 



Levy fc Lunind Il993l . iMannings. Boss &: Russel] 2000l . iReipur th. Jewitt k. Keil 



20071 ). There are also several related reviews that have been rec ently written 



for this series including the inner disk toullemond &: M onnier 20ld ) , debris disks 



( Wvatt|[2008l ). and dynamical processes (|ArmitagQj201Q l 



This review focuses on the properties and evolution of the outer parts of pro- 
toplanetary disks around low mass stars as determined principally from obser- 
vations at mid-infrared to millimeter wavelengths. After briefly describing the 
classification of young stellar objects in ^ we begin by discussing the formation 
of disks in ^ and their basic properties when the central star first becomes opti- 
cally revealed in ^ The second half concerns itself with the temporal properties 
of disks. We discuss their lifetimes in ^ the evidence for, and processes by 
which, disks evolve in ^ and the properties of disks transitioning into their end 
state in ^ Each section ends with a short summary of the key points and these, 
in turn, are distilled into an overall summary in ^ There are many promising 
avenues for future exploration that new and planned facilities can address, and 
we list those that we deem most exciting in ^ 

2 CLASSIFICATION OF YOUNG STELLAR OBJECTS 

The process of star and planet formation begins with the collapse of a molecular 
core. The mass is initially all in the core but it is processed through an accretion 
disk inward onto the protostar and outward through an outflow. Ultimately the 
core is dispersed and the remaining mass is concentrated in the star. There 
i s enough terminolo gy associated with this process to form its own "diskionary" 
( Evans et al ] l2009al ) and many observational ways to characterize the progression. 



The most direct, that of measuring the mass in each component, is hard and 
more practical means are generally used, principally measuring disk accretion 
signatures in the optical or the distribution of warm circumstellar material in the 
infrared. 

The infrared based classification dates back to Lada k. Wilking ( 19841 ) who 



showed that Young Stellar Objects (YSO) in Ophiuchus formed 3 distinct groups 
based on whether the emitted energy was rising in the mid-infrared, declining but 
with a notable excess over the blackbody stellar photosphere, or with negligible 



infrared excess. This was formalized into 3 classes, I-II-III respectively, by iLada 
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(|l987l l based on the slope of the spectral energy distribution (SED) between 
about 2 and 25 /um, 

dloguF^ dlogXFx 



am 



d log u d log A 

Greene et all (|l994l ) subsequently introduced an additional refinement of "flat- 



spectrum sources", intermediate between Class I and II YSOs. The class sequence 
was sho wn to fit naturally into the th eoretical framework of a rotating, collapsing 
core by Adams. Lada Shu ( 19871 ). As the ability to detect faint millimeter 
em ission improved, the categorization was ext ended to an earlier, Class 0, phase 
by Andre. Ward-Thompson &: Barsony ( 19931 ) and the decrease in mass of the 
circumstellar envelope (m aterial participating in the star/disk formation) along 
the sequence was verified ( Andre fc Montmerlel[l993 ). 

A parallel accretion-based classification exists for the later optically visible 
phases: classical and weak-lined T Tauri stars. These correspond closely (though 
not exactly) to Class II and III YSOs respectively. Classical T Tauri stars (CTTS) 
have strong Ha and UV emission whereas weak-lined T Tauri stars (WTTS) show 
no or only very low indications of accretion. Historically the dividing line between 
the two was a uniform Ha equivalent width of lOA but this has been refined to a 
stellar mass dependent limit to account for t he lower continuum l evel in low mass 
stars ( Barrado v Navascues &: Martini 20031 . I White k Basrill2003l ). The census of 
Class HI sources or WTTS is generally incomplete in optical and infrared cata- 
logs because they lack strong accretion signatures and infrared excesses. Their 
youth is generally only betrayed through their location in the Hertzsprung-Russell 
diagr am above the main sequence or by X-ray activity (jFeigelson &: Montmerle 
19991 ). 



It is also important to note that the SED classification does not give a unique 
description of the amount and distribution of circumstellar material. In particu- 
lar, YSOs with edge-on disks are highly extincted and can be mis-interpreted as 
more embedded, hence, less evolved, objects. For example, a Class II YSO viewed 
at high inclination has a similar SED to a typical Class I YSO and an edge-q n 



Class I YSO can have characteristics of a Class YSO (iRobitaille et al.ll200(il ) 



These ambiguities highlight the necessity to distinguish between the physical 
"stage" of YSO evolution for comparison to theory and its observed SED. Re- 
solved images, ideally at multiple wavelengths, are required to fully characterize 
the evolutionary state of any individual YSO. Table 1 sumn iarizes the definition s 
for each class, where the numerical boundaries for aiR follow Greene et al. ( 19941 ). 
including corresponding physical properties (evolutionary stage) and other obser- 
vational characteristics. 

A disk forms very early on and grows rapidly during the Class collapse phase 
(see jQl). On average, the e mbedded phases through Class I lasts for about ~ 
0.5 Myr ( Evans et al.l 2009bl ). The surrounding envelope and strong protostellar 
outflows can hinder measurements of disk properties at these early times. The 
properties of the disks as their central stars first become optically visible, i.e. 
Class II YSOs, are discussed in ^ The median disk lifetime after the embedded 
phase is 2-3 Myr but the manner and the rate at which any individual star-disk 
system evolves vary greatly. These issues are discussed extensively in §^[6l 

Large Spitzer surveys h ave mapped about 90% of all the star- forming regions 
within 500 pc of the Sun ([Evans et al.ll2009bl . and references therein) and spec- 



tra have been obtained for over 2000 YSOs therein (e.g.. iFurlan et al.l I200S 



Kessler-Silacci et al.l l2007l . lOliveira et aP l20ld '). The precision and wavelength 
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coverage of these observations show the tremendous diversity of disk SEDs and 
the inadequacy of a single parameter, ajR, to characterize the full range, espe- 
cially as disks dissipate and open up central holes. These so-called transition 
disks are the subject of ^ 



3 DISK FORMATION 



The initial collapse of a molecular cloud core is onto a point source but a disk 
quickly forms as more distant material with higher angular momentum falls in- 
ward. The disk extends out to the centrifugal radius, which is expected to grow 
rapidly with time, R(t) oc Q'^t ^, where Q is the angular rotation rate of the core 
( Terebev. Shu &: Cassen 19841 ) . Disks should evolve rapidly, therefore, and their 
final size and niass dep ends sensitively on the infall time (t^) and the core proper- 
ties (ri^). iBasu ( 19981 ) notes that magnetized collapsing cores may not be in rigid 
rotation and that the radius may grow only linearly with tirn e. In any case, given 
the wide range of core rotation rates ( Goodman et al. 19931 ) and likely variation 
in infall duration, we should not be surprised by an inherent and large diversity 
in initial disk sizes and masses. 

The role of magnetic fields in core collapse and disk formation is uncertain. Po- 
larization measurements show magnetic field lines concentrat ing at the center of 
collap sing cloud cores in a pinched, "hour glass" configuration (jGirart. Rao &: Marrone 
[2006') yet Zeeman observations show that the magnet ic field strength is generally 



insuffi cient to support cores against their own gravity (ICrutcher. Hakobian &: Troland 

We expect some field lines to be dragged down to disk scales for the 
magneto-rotational instability, the most likely mechanism for disk viscosity at 
late ti mes, to exist. Sing ledish sub-millimeter measurements of polarization in 
disks (jTamura et al.lll999l i have not been co nfirmed with interfero metry, however, 
thus challenging theoretical understanding ( Hughes et al. 2009bl ). 

Numerical models of both magnetic and non-magnetic collapsing molecular 

cores show disks form rapidly, within ~ 10^ yr ( Hueso &: Guillo'OboOSi : Yorke. Bodenheimer &: Laughlin 
19931 ). Temperatures are very high in these early stages due to gravitational in- 



fall. As the core material is used up or otherwise dispersed, the disk cools down 
and its mass decreases as it accretes onto the star. 

Core collapse onto a disk opens up an approximately spherical cavity in the 
surrounding envelope of radius R{t) that has been inferred from the presence 
of excess mid-infr ared emission above that expected from a more extincted cen- 
trally peaked core ( Enoch et al.ll2OO9l . |j0rgensen et al. Il2005b|). Although there are 
many observations of inward motions on core-size scales (e.g.. lDi Francesco et al. 
200ll ). the dire ct detection of g as flo w onto a disk has yet to be convincingly 



demonstrated. Chandler et al 



(|2005) find SO line absorption against the disk 
continuum in the IRAS 16293-2422B YSO but the spectral profile is approxi- 
mately sym metric about the sou rce velocity and cannot be clearly identified as 
pure infall. IWatson et all (l2007l l detect many mid-infrared lines of H2O toward 
NGC1333-IRAS4B, which they model as arising from a dense, warm, and compact 
region. T hey attribute this to shocked gas from an envelope onto a disk surface. 
However, Ijgirgensen van Dishoeckl (|20inl ) mapped the 1.5 mm 3i,3 — 22,0 transi- 
tion of H2^0 from the ground and find the water emission is quiescent and follows 
the disk rotation. 

Imaging embedded disks requires long wavelengths to see through the en- 
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velopes and arcsecond or higher resolution to match the disk sizes. Milhme- 
ter inteferometers meet these requirements and also filter out extended struc- 
tures so that the emission from the compact disk domi i iates on long baselines , 
> 50kA (e.g ■ iBrown et al."2000l. Ijoreensen et alJl2005al . iKeene Massonlll990l . 
Looney. Mundv &: Welch 2000). A 1.1mm continuum survey of 20 embedded 



YSOs bv Ijgreensen et aD (|2009) shows that the disk flux is an average of four 
times higher in the embedded Class phase than in Class I sources but, after 
allowing for higher temperatures due to greater accretion heating, the inferred 
disk masses show no significant dependence on evolutionary state. Masses in both 
Class and I sources range from ~ 0.02 — 0.1 Mq with a median 0.04 M©. This 
is quantitatively similar to the results from 2.7 mm observations of 6 sources by 
Loonev. Mundv k WelchI (|2003l l. 



The lack of dependence of disk mass on evolutionary state from Class to 
Class I is contrary to theoretical expectations of steady disk growth as outlined 
above, whether the core is rigidly or differentially rotating. Rather, it indicates 
that disks form quickly and that the flow of material from the envelope, which 
declines in inass b y almost an order of magnitude between these two classes 
(jYoung et al.ll2003l ). is rapidly transported through the disk. 

A likely cause of the rapid tr ansport at these early stages is disk instability. 
Laughlin &: Bodenheimer (| 19941 1 first suggested that disks would be gravitation- 



ally unstable during the early formation stages due to the relatively high mass 
fraction in the disk versus that accreted onto th e protostar. Th e instabilities 
would lead to sporadic bursts of high accretion (IZhu et al.l l2009l'l and prevent 
the disk mass from growing faster than the star ( Vorobyov &: Basu 2O10l ). Many 
protostars have been observed to undergo short-lived bursts of activity due to 
high acc retion and th es e events are named after the first such identified case, FU 
Orionis (lHerbigi[l977l V iHartmann fc KenvonI (|l99()l ) estimate that a typical low 
mass star may have an average of about 10 such outbursts during its formation. 

Additional supporting observational evidence for "punctuated evolution" at 
these early stages is found in measurements of mass infall rates through the 



differe nt components and the protostellar luminosity distribution. lEisner et al. 
(I2OO5I ) find envelope infall rates are more than an order of magnitude higher than 
disk accretion rates in Class I YSOs suggesting that mass builds up in the disk 
until a burst event occurs. It has long been known that Class I YSOs are, on 
average, about an order of magnitude less luminous than expected for the steady 
release of gravitational energy as the envelope fall s onto the protostar over the 
lifetime of the embedded phase ( Kenvon et al. 199d ). Based on the st atistics from 



the cor es-to-disks Spitzer survey of 5 large, nearby molecular clouds, lEvans et al. 



(12009^ ) concludes that, on average, a star gains half of its final mass in only 
7% of the ~ 0.5 Myr Class plus I lifetime. These statistics c an be matched 



by m odels of episodic mass accretion from self-gravitating disks (jPunham et al. 
2OIOI I. 



Some words of caution are necessary in the interpretation of the interferometric 
data, however. The analyses of the continuum visibilities are based on rather 
sparse sampling of the Fo urier plane and are s ubject to confusion with small scale 
structure in the envelope (IChiang et al.ll2008l ). A more secure identification of an 
embedded disk requires spectral line observations showing rotation. The emission 
from the core and any protostellar outflow may still present significant confusion, 
although this is somewhat mitigated by observing species such as HCN and HCO"*" 
that only emit substantially in gas that is warmer and denser than the outer parts 
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of the envelope. Recent work demonstrates the possibihties in this area with 
i ndications of Keplerian velocity profiles in moderately young , Class I sources 



( Brinch et al. 2007 . Jorgensen et al.l 20091 . Lommen et al. 20081 ) . More sensitive 



observations of optically thin isotopologues are required to discern disk kinematics 
from the core background in the younger, more embedded Class phase. Such 
observations also have the potential to measure the central protostellar mass and 
track its growth against the evolution of the envelope and disk. 



3.1 Section summary 

• Circumstellar disks form almost immediately after a molecular core col- 
lapses. Their existence can be inferred through both SED modeling and 
interferometry. 

• Disk masses do not appear to increase with time during core collapse im- 
plying a rapid transport onto the star. 

• The low average luminosity of YSOs and the accretion bursts of FU Orionis 
objects suggest that young disks are (gravitationally) unstable. 



4 PROPERTIES OF PROTOPLANETARY DISKS 



The deeply embedded Class and I phases of star formation only last for a 
small fraction of a disk lifetime, typically ~ 0.5 Myr compared to several Myr. 
By the end of the Class I phase the envelope has completely dispersed and the 
star formation process is effectively over. The disk now contains only a few 
percent of the central stellar mass and can be considered truly protoplanetary, 
not protostellar. Although there may be a small amount of accretion of material 
from the molecular cloud, the major elements governing the evolution of the disk 
at this stage are accretion onto the star, photo-evaporation from local or external 
radiation sources, agglomeration into larger bodies, and dynamical interactions 
with stellar or substellar companions. 

Unless the disk is edge-on to our line of sight, the extinction to the central 
protostar is small and its spectral type can be determined through optical or 
near-infrared spectroscopy in most cases. The properties of the disks at this 
stage marks the baseline for studies of their evolution and for core accretion 
models of giant planet formation. In this section, we discuss the basic properties 
of the outer regions of protoplanetary disks around Class II YSO. 

As a fiducial comparison, we use the Minimum Mass Solar Nebula (MMSN), the 
lowest mass primordial disk that formed the Solar System inferre d from scaling 



plane t ary com position s to co smic abundances at each orbital radius (jKusaka. Nakano &: Hayashi 
I970I : Weidenschilling 1977). Becau s e of t he uncertainty in the composition 



of the giant planets, Weidenschilling ( 19771 ) quotes a range for the MMSN of 
0.01 — 0.07 Mq. Here, we use the minimum of this range, 0.01 Mq w lOMjup, as 
an absolute lower mass for the solar nebula out to the 30 AU orbit of Neptune. 
Extrasolar planetary systems have been found around 10% of Sun-like stars in 
the solar neighborhoo d, and m any planets detected to date are considerably more 
massive than Jupiter ([Johnson 2009). The minimum mass disk required for their 
formation is proportionally higher. 

The extrapolated surface density profile, which is an observable for resolved 
disks, has an approximate power law form, S oc r~^/^, although more sophisti- 
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cated fits based on viscous dis k evolution and planet migration models can be 
made (|Davidl2005l . lDeschll2007l l. 

The properties of the outer parts ( ^ 1 AU) of protoplanetary disks are mainly 
inferred from observations at mid-infrared through millimeter wavelengths al- 
though optical and near-infrared scattered light or silhouette obser vations can 



provide important inforni ation on disk radii and dust grain properties (jPinte et al, 
20081 . iThroop et alJ[2nnil ). 



4.1 Mass 

Disk masses are best determined from (sub-)millimeter wavelength observations 
of dust. The continuum emission is optically thin except in the innermost regions 
where column densities are very high. The optical depth is the integral of the 
dust opacity, k^, times the density, p, along the line of sight, Ty = ^ pnyds = k^Yj, 
where S is the projected surface density. A commonly used prescription for the 
dust opacity in disks at millimeter wavelengths is 



0.1 



1012 Hz 



2 -1 
cm g 



(2) 



(|Beckwith et al.lll990l ^. Both the absolute value and power la w index, /3, are re- 
lated to the size distribut ion and composition of the dust grains ( Ossenkopf &: Henningj 
1994 , Pollack et al.|[l993 ) . The normalization above also implicitly includes a gas- 
to-dust ratio of 100, and p and S refer to the total (gas plus dust) density. 

For disks, /3 w 1 so k(1 mm) = 0.03 cm^ g ^ which implies r(l mm) = 1 at a 
surface dens ity of S w 30 gcm"^^ corresponding to about 10 AU in the MMSN 
([Davisl I2OO5I ) and an angular scale of 0. 07 in the nearby Taurus star-forming 
region. Protoplanetary disks are generally much larger than this and most of 
the resolved emission is indeed optically thin. However, the inner 10 AU likely 
constitutes a substantial fraction of the planet-forming region of a disk. Not only 
high resolution but also wavelengths longer than 1 mm are required to peer into 
this zone at the early phases of disk evolution. 

For the general question of disk mass measurements on scales much larger than 
10 AU, we can consider the emission to be optically thin and therefore directly 
relate the observed flux, Fy, to the mass, 



M(gas + dust) 



(3) 



where d is the distance to the source. This shows an additional advantage of 
millimeter wavelengths in that the Planck function is close to the Rayleigh- Jeans 
regime. By « 2v^kTl(?, and the emission is only linearly, rather than exponen- 
tially, dependent on the dust temperature. 

The first large millimeter wavelength surveys to measu re disk masses w e re car - 
ried out at 1.3 mm usi ng a single element bolome ter by Beckwith et al. ( 1990l ). 



in Taurus- Auriga and Andre &: Montmerld ( 19941 ) in Ophiuchus. Both regions 
are predominantly forming low mass stars, M^, < 1 Mq with spectral types K-M. 
Using bolometer arrays, which provide better sky subtrac tion capability, hence 
l ower noise l evels and the ability to detect fainter targets, [Andrews &: Williams 
(|200,1 l2007al ;i augmented these surveys in size and expanded the wavelength cov- 
erage into the sub-millimeter regime with observations from 350 pm to 850 pm. 
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The results of this work show that the median mass of Class II YSO disks is 
5Mjup and the median ratio of disk to stellar mass is 0.9%. Figure [1] shows 
the mass distributions are rather flat by logarithmic mass interval with a steep 
decline beyond ~ 50Mjup. The e xplicit labeling of th e gas-to-dust ratio in this 
figure mimics an earlier review by lBeckwith fc SargentI (Il996i ) and serves to illus- 
trate the large extrapolation necessary to derive total disk masses and constrain 
models of giant p l anet f ormation. This figure also includes the disk census by 
Mann fc Williamsl ^201(h in the more distant but massive star forming Orion 

region where photoevaporation is important fsee §5. 6. 2D. 

By modeling the infrared-millimeter SED, Andrews &: Williamsl ( 20051 ) showed 
that the simple mass derivation in equation [3] is indeed reasonably accurate with 
a characterist ic temperature , T = 20K. This temperature i s consistent with CO 



obser v ations (iQi et al.l 120041') and t heoretical expectations (jChiang Goldreich 



19971 ). [Andrews fc Williamsl (|2005l ) also show that for a 100 AU disk with the 
mass and surface density profile of the MMSN, about one third of the emission 
by mass is optically thick at 850 /im. 

An additional major uncertainty is the hidden mass in large grains. Dust grows 
to much larger sizes in protoplanetary disks than in the interstellar medium and 
these "pebbles" or "snowballs" can hold considerable mass in a small solid an- 
gle with negligible effect on the SED. As a rule of thumb, observations at a 
wavelength A only constrain the properties of dust grains out to a maximum 
size Cmax ~ 3A (Draine 2006). Only a few disks have only been detected be- 
yond millimeter wavelengths, and thus we know little about the general occur- 
rence and distribution of centimet er and larger sized particles in dis ks. For a 
grain size distribution n{a) oc a~^'^ ( Mathis. Rumpl &: Nordsieck 19771 ). the total 

1/2 

m ass scales as Omax and substantia l mass may be undetected. Detailed modeling 
by D'Alessio. Calvet &: Hartmann (|2001i ) shows that the opacity prescription in 
equation [2] is valid for Omax ~ 0.3 — 3 mm but is about a factor of 20 smaller for 
Omax = 1 m implying correspondingly higher masses. The observational evidence 
for grain growth from sub-micron to millimeter sizes and beyond is discussed in 

The uncertainties in disk mass measurements are large but fortuitously tend to 
cancel each other out. We overestimate by assuming an interstellar gas-to-dust 
ratio and underestimate by ignoring large bodies. There are a couple of indica- 
tions that disk masses as derived from equations [2| and |3| are significantly low, 
however, which would imply that grain growth is the dominant effect. First, pro- 
toplanetary disk masses derived from (sub-)millimeter photometry are systemati- 
cally lower than those estiinated f r om accretion rates inte grated over protostellar 
ages (jAndrews &: Wilhamsl 12007^ . iHartmann et al.lll998l ). Second, there are not 
enough massive disks in nearby star-for niing regions to match the statistics on 
the incidence of extrasolar giant planets ( Greaves &: Ricell20ld ). 



4.2 Radius 

Disk sizes are hard to measure because the outer parts are cool and emit weakly. 
They are still efficient absorbers, however, and images of disk silhouettes in Orion 
against the opti cally bright H II region background provide a simple and direct size 
determination. IVicente &: A lves (200^) measure radii ranging from 50 to 194 AU 
for 22 Orion "proplyds" where the disk shadow can be clearly seen in Hubble opti- 
cal images and an additional two outliers with radii of 338 and 621 AU. The latter 
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object, cataloged as proplyd 114-426, is almost 3" in angular extent, resolved in 
both radial and vertical dimensions, and provides perhaps the most famous dis k 
image obtained to date ( McCaughrean et aPllQQSl . lMcCaughrean fc 0'dell|[l99^ ). 
Most of the Orion proplyds are ap parent only throu g h the photoevaporative flows 
from their surfaces, however, and Vicente &: Alves ( 20051 ) estimate radii for 125 
such objects from the size of the ionization front. Although there is substantial 
uncertainty associated in this indirect determination, they infer a median radius 
~ 75 AU. They further note that the sample of Hubble proplyds accounts for only 
half the stars with known infrared excesses and suggest that over three quarters 
of the full sample of Orion disks have radii less than 75 AU. 

Imaging disks at millimeter wavelengths requires intefer ometry on acc o unt o f 



their small angular scales in nearby star-forming regions. iDutrev et al.l (|l996l ) 
carried out the first large interferometric survey and resolved many Taurus disks 
with typical angular sizes of 1 — 2". A curious problem arose in that the size of the 
gas disk, as observed in rotational lines of CO, was found to significantly exceed 
the size of the continuum image (jlsella et al.ll2007l : IPietu. Guilloteau fc Dutrev I 
20051 ). The differences could not be reconciled for a sharply truncated power law 
surface density profile (i.e., S oc for R < Rout, S = for > Rout) without 
introducing an arbitrary change in either the d ust-to-gas ratio or dust o pacity 
at the radius of the continuum disk. Similarly, McCaughrean &: O'delll ( 19961 ) 



found that pure power laws or sharp edges do not match the intensity profiles 
of Orion proplyd silhouettes and showed that an exponential decay at the outer 
boundary was required. Indeed physical models of viscous accretion disks (e.g., 
Hartmann et al.l 19981 . Lynden-Bell &: Pringle 19741 ) , predict an exponentially ta- 
pered profile of the form. 



exp 



R 

Rr 



2-i 



(4) 



where is the disk mass, Rc is a characte ristic radius, and 7 spec ifies the radial 
dependence of the disk viscosity, v oc IV . iKitamura et al. I (|2002l ) first modeled 
millimeter wavelength disk images with this profile, but their data were unable 
to differentiate between this and sharply tru ncated power law fits. With higher 
resolution observations, Hughes et al.l ( 2008 ) showed that this prescription can 
naturally account for the apparent size discrepancy in millimeter imaging as the 
the optically thick CO line can be detected further out along the tapered edge 
than the optically thin dust continuum. 

Rc is a characteristic radius that delineates where the surface density profile 
begins to steepen significantly from a power law. As the disk does not have a sharp 
edge, a physical disk size must be specified in terms of an intensity threshold. A 
rough estimate of Rc may be obtained by noting that abou t 2/3 of the total disk 



mass (approximately flux) lies within it. More precisely, Hughes et al. (j2008l ) 
determine i?c = 30 — 200 AU for four disks by fitting the profile in equation |4] 
to the continuum data. They compare with sharply truncated power law fits 
and find Rout ~ '^Rc- The tapered disks extend well beyond Rout, however, with 
midplane densities n(H2) > 10^ cm~^ for R ^ 800 AU. 

By s imultaneously fi t ting 850 jum i nterferometric visibilities and infrared-millimeter 
SEDs, [Andrews et aD (l2009l . l2010bl ) determine a wide range Rc = U - 198 AU 
for 16 disks in Ophiuchus. Their sample is sufficiently large that they are able 
to identify a correlation with disk mass, oc R^'^^^'^, which is likely to be a 
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signature of disk formation physics rather than an evolutionary sequence. They 
d o not find any correlation wi t h the central star properties. 

Isella. Carpenter &: Sargent ( 20091 ) use the same density prescription to fit a 



sample of 11 disks mostly in Taurus- Auriga. They use a slightly different radial 
normalization that converts to a similar range, Rc — 3 — 230 AU, and also 
shows that the older disks tend to be larger. In addition, Schaefer et al. (j2009l ) 
compile radii for 16 disks around low mass stars in Taurus- Auriga based on 
sharply truncated power law fits to CO observations. Because of the different 
fitting technique, their values. Rout ~ 100 — 1100 AU, are not directly comparable 
with the above Rc measures, but they similarly illustrate the large range of disk 
sizes and independence on stellar type. 

The specific angular momenta o f protoplanetary disks, j = J/M, lies in the 

range logj^Q j [cm s" -1] = 19.4-20.9 (|Andrews et al.ll2010bl : llsella. Carpenter Sargent! 
20091 ). This is comparable to the giant planets in the Solar System, log;^oJ = 
20.0 — 20.4, but signi ficantly lower than th e observed values in molecular cores, 
logio j ~ 20-^ ~ 22.6 ( Goodman et al. 19931 ). Either the angular momentum in a 



core is not efficiently transferred to the disk, or its value is overestimated due to 
the inherent sin oothing of turbulent fluctuations in radial velocity observations 
(|Dib et al.ll2010l l. 



4.3 Structure 

The average radial, vertical, and velocity profiles of relatively bright and large 
disks can be determined through resolved observations and careful modeling of 
spectral energy distributions. 

4.3.1 Surface density A resolved image of a disk at millimeter wave- 
lengths provides not only a measure of its total mass and radius but also the 
distribution of mass, or surface density. Until recently this was characterized 
as a pure power law, S oc R~^, with values of p generally in the range — 1 
(lAndrews fc Winiamsll2007bl. iKitamura et al.l [2002. iLav. Carlstrom fc Hillsl[l997l . 
Mundv et al.lll996l . IWilner et al.ll20o'nf K 

The exponential tapered fits of the form in equation |4] approximate a power 
law, S oc R~'^ for R <C Rc, but the fitted values of i?c — 30 — 200 AU correspond 
to ^ 1" in all but the closest disks and 7 is actually determi ned largely from the 
steepness of the exponential taper. The Hughes et al. ( 20081 ) comparison of pure 
power law versus exponentially truncated power law fits shows similar indices but 
with slightly steeper pure power law fits due to the soft edge, (p) = 1.2, (7) = 0.9 
for four disks. 

From larger samples, Andrews et al. (2009, 2010bl ) find a tight range consistent 
with all having the same value (7) = 0.9. U sing a different modeling technique, 
however, Isella. Carpenter &: Sargent ( 20091 ) find a very wide range 7 = —0.8 
to 0.8 with mean (7) = 0.1 in their data. Negative values of 7 correspond to 
decreasing surface densities for R < Rc which may be an important signature 
of disk evolution (see ^ but these were also found in the smaller disks with 
Rc < 100 AU which are barely resolved and thus the hardest to characterize. 

In general, all results agree that young protoplanetary disks ha ve flatter central 
densi ty proflles than the canonical power law p = 1.5 MMSN ( Weidenschillingj 
19771 ). There is even more uncertainty in the density proflle of the MMSN than 
its mass, however, and an exponential tapered power law fit bv lDavisI (|2005l ) has 
7 = 0.5. The more relevant comparison is of absolute values in the planet-forming 
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zone. Andrews et al. I (|2009l . l2010bl 'l infer surface densities, S ~ 10 — 100 gem 
at 20 AU, in their sample of Ophiuchus disks that are in good agreement with 
the MMSN (Figure [2]). 

Whereas the disk-to-star mass ratio may be very high during the initial stages of 
formation the Toomre Q paramete r, Q(R) = cQ/ irGTi where c is the sound 
speed and J7 the orbital angular velocity ( Toomrg 19641 ) . is generally much greater 
than unity for th ese Class II YSO disk s implying that they are grayitati onally 
stable at ah radii ( Andrews et al.|[2010bl : llsella. Carpenter Sargentll2o"09l ). The 
one possible exception is DC Tau, a young protostar with by far the most massive 
disk in the two surveys. 

The inferred surface densities are increasingly uncertain closer to the star due 
to the limited resolution of the observations, ^20AU, and al so because the 



emission is be comi ng optically thick. As the discrepan cy between [Andrews et al. 



(|2009l . l2010bl ) and llsella. Carpenter k SargentI (|2009l ') shows, different data and 
fits extrapolate to very different surface densities at the 5 — 10 AU orbital radii 
of Jupiter and Saturn. This is not only a critical region for understanding planet 
formation but also for disk structure. The density is so high here that the mid- 
plane is shielded from radiation and cosmic r ays. Unless radioactive elements are 
sufficiently abundant ( Turner fc Drakj 20091 ). the ionization fraction may be so 
low that matter is decoupled from the magnetic field and avoids the magneto- 
rotational instability (see ^6.ip . Current observations cannot address these issues, 
and they will remain a challenge until the combined requ irements of h igh sen- 
sitivity and high resolution at long wavelengths are met (IWilneil[2003 ). Note, 
however, that the uncertainty in the surface density is outweighed by the Kep- 
lerian increase in rotation rate and Class II YSO disks are gravitationally stable 
in their inner regions. 

4.3.2 Scale height Proto planetary disks are flared w ith a vertical scale 
height that increases with radius. iKenyon &: Hartmann (Il987l ) first suggested the 
possibility of flaring based on the large far-infrared excesses detected by IRAS 
that could not be explained by a spatially flat disk. Direct evidence for flared disks 
can be found in the beautiful Hubble images of Taurus disk silhouettes against 



the scattered light of the central star (jBurrows et al.lll996l . iPadgett et al.l 11999 



Stapelfeldt et al. 19981 ) and in Orion against the nebular background ( Smith et al. 
20051 ). 



Characterizing the disk scale height is essential for modeling the thermal, ion- 
ization and chemical structure of disks and, thence, for interpreting atomic and 
molecular line observations. It is also important for understanding disk evolution 
as the outer parts, with their low densities and large scale heights, are particularly 
susceptible to photoevaporative losses. 

For an azimuthally symmetric disk in hydrostatic equilibrium, the density is a 
function of both radius, R, and vertical height, Z, 



p{R,Z) 



2ttH 



exp 



2IP 



(5) 



where S(i?) is defined in equation IH and H{R) is the scale height which depends 
on the competition between disk thermal pressure and the vertical component of 
stellar gravity. The disk temperature, in turn, depends on the amount of stellar 



radiat ion impacting the disk and therefore on its geometry. IChiang &: Goldreich 



(119971 ) provide an elegant analytic solution to these coupled equations by ap- 
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proximating the radiative transfer with a hot surface layer that absorbs and 
reprocesses the starhght, which then heats a flared interior. T hey find an ap- 



D'Alessio et al. 



proxim ate power law dependence, H oc R^, w i th h 1.3 — 1.5 

and lDullemond. van Zadelhoff fc Nattal \2m± iterate numerical solutions 
to the vertical disk structure and find similar results. 

Observed disk SEDs have less mid-infrared emission than exp ected for this 



degree of flaring due to sett ling of dust grains toward the midplane (IChiang et al 



2nnil . iD'Alessio etld]ll999l ). This flattening of the dust disk relative to the gas 



occurs very early in the evolution of a disk and is discussed further in ^ 

4.3.3 Velocity Disk masses at the Class II stage and beyond are a small 
fraction of the central stellar mass. Their motions are therefore expected to be 
Keplerian. Velocity profiles have been measured for only a handful of disks. The 
difficulty lies in finding disks that are sufficiently young and therefore bright 
enough to image in a spectral line, generally a millimeter wavelength rotational 
line of CO or isotopologue, yet not contaminated by emission from the residual 
ei ivelope or neighboring cloud. 

Koerner. Sargent &: Beckwith ( 19931 ) resolved the ^^CO J = 2 — 1 line in the 



relatively large CM Aur di sk and found the velocity profile w as consistent with 
being Keplerian. Similarly, iDutrey. Guilloteau &: Simon imaged the large 

circumbinary disk around GG Tau in ^^CO J = 1 — and used a Keplerian 
velocity profile to determine a total stellar mass of 1.2 Mq. With improve- 
ments in instrumentation, leading to higher frequency an d higher resolution 



1998; 



observations, more studies followed shortly thereafter (e.g., iDuvert et al. 
Mannings. Koerner fc Sargentlll997l ). 

Guilloteau k. Dutrey (jl998l ) modeled the CO J = 1 — channel maps (images 



at different velocities) of DM Tau to show that a Keplerian velocity profile was 
no t only consistent but was the b est fit . This was further exploited in the survey 
by Simon. Dutrey &: Guilloteaul ( 2000l ) who determined dynamical masses for 9 
systems and tested models of protostellar evolution. Schaefer et al. (2009) ex- 
tended this work to later spectral types. As discussed in ^ Lommen and Brinch 
have used the same technique to derive protostellar masses at even earlier times 
in Class and I YSOs. Spectral line observati ons also determ ine disk inclinations 
to our line of sight with high accuracy (e.g., Qi et al.l 20041 ) . which is necessary 
for detailed modeling and characterization of other disk properties. 

Many of the aforementioned studies also note that the line emission is well fit 
by only rotational plus thermal broadenin g. High spectral reso lution observations 
of the TW Hydra and HD 163296 disks bv lHuehes etall (|2010l ) show that the tur- 
bulent component is subsonic, < 10% and 40%, of the sound speed respectively 
at the ~ 100 AU scales of their observations. Such a low level of turbulent stirring 
provides ideal conditions for grain settling and growth, planetesimal congrega- 
tion, and protoplanetary accretion. Their data also show near-perfect Keplerian 
rotational profiles confirming that disk self-gravity does not appear to be signifi- 
cant in these objects (Figure [3]). 

4.3.4 Departures from azimuthal symmetry The discussion in this 
section has so far considered only the radial or vertical variation of disk properties 
on account of the dominant gravity of the central star. Azimuthal variations are 
of great interest, however, because they would signpost an additional influence, 
potentially disk self-gravity or proto-planets. 

Disks are tiny compared to the distances between stars even in dense clusters 
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and disruptive encounters are highly unlikely. For typical relative motions of 
1km s~^, the stellar density must be greater than ~ 10^ stars pc~^ for an en- 
counter within 1000 AU of a disk within 1 Myr. Almost all disks that have been 
studied in detail are nearby and are either isolated or in low mass star-forming 
regions where the stellar densities are orders of magnitude lower than this. Con- 
sequently, the likelihood that any azimuthal distortions in a young disk are due 
to an encounter with a nearby star is very low. 

The surface density measurements described in §4.3.11 and Keplerian veloc- 
ity profiles in ^4.3.31 show that disk self-gravity is negligible by the time of 
the Class II phase. As discussed in ^ however, younger disks in formation 
probably have much higher disk-to-star mass ratios. The disk around the par- 
tially embedded Herbig Ae star AB Aur may be an exai nple of an u n stable 
disk . Optical and ii e ar-inf rared coronagraphic imaging bv iGradv etld] (|l999l l 
and Fukagawa et al. ( 20041 ) show spiral features in scattered light. Lower res- 
olution millimeter interferometry confir ms the larger scale asymmetries in the 
system ( Corder. Eisner &: Sargent 20051 ) and deviations from Keplerian rotation 
that may indicate streaming motions along the a rms or an unseen companion 
(|Lin et al.ll2006l : IPietu. Guilloteau Dut"revll2005l ). 

The gravitational effect of planets on debris dis ks has been diagnose d in Foma- 



Ihaut (jChiang et al.1 120091'). p ostulated in /3 Pic (jMouillet et all 119971 ). and sus- 



pected in others ( Wvatt 20081 . and references therein). Younger disks are further 
away than these systems and, the intriguing case of AB Aur notwithstanding, 
the resolution is not yet sufficie nt to search for such effects, although this is a n 
exciting prospect for the future ( Narayanan et al.ll2006l . IWolf &: D'Angeloll2005l ). 



4.4 Composition 



As with the interstellar medium, the circumstellar medium is composed of dust 
and gas. There is substantial processing of both components, however, as a disk 
forms and evolves. 

4.4.1 Dust Dust dominates the opacity of protoplanetary disks and is the 
raw material for planetesimals. In the diffuse interstellar medium, dust is mainly 
composed of silicates with sizes r ^ 0. 1 um with an admixture of graphite grains 
and polycyclic aromatic hydrocarbons (|Drainell2003l V During the passage through 



cold dark clouds to protoplanetary disks, molecu 
onto dust grain surfaces produci ng icy mantles 



grains collisionally agglomerate (jBlum fc Wurm 



es freeze out from the 



Bergin k Tafallall200 



20081 ). 



;as phase 
and the 



Silicates are readily identified through broad spectral bands at 10 and 18 ^m 
( Henninj 2010l ). Polycyclic aromatic hydrocarbon s also have several mid-infrared 
features in this wavelength range (|Tielensl[2008l l. About 2000 protoplanetary 
disks were observed with the Spitzer Infrared Spectrometer (IRS) providing a 
wealth of information on the composition, size distribution, and evolution of the 
dust. We defer discussion to ^6.31 and note here only that the dust in disks has 
been significantly processed via thermal annealing and agglomeration to have 
a much higher crystallinity fraction and larger sizes relative to the interstellar 
medium. 

4.4.2 Gas Gas amounts to 99% of the total mass of the interstellar medium 
and the same is true, at least initially, for disks. Its detection, however, presents 
a tougher observational challenge than the minority dust component because it 
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emits only at specific wavelengths and therefore requires high resolution spec- 
troscopy. Moreover, at wavelengths where the dust emission is optically thick, a 
line can only be seen if its excitation temperature differs from that of the dust. 

Disk accretion, which can be measured through recombination line emission or 
excess hot continuum emission, provides unambiguous evidence for the presence 
of gas. However, the precise amount or physic al conditions of t he gas in the disk 
cannot be determined from these diagnostics ( Hartmann 20091 ). 

At the typical high densities and low temperatures of the bulk of a protoplane- 
tary disk, the primary gaseous species is H2. Fluorescent el ectronic transitions o f 
H2 in the UV arise from hot gas very close to the protostar ( Herczeg et al. 20061 ). 
Near-infrared ro- vibrational lines also probe the inner disk, ^ 1 AU. Because the 
focus of this review is on the o uter disk, we d iscuss mid-infrared observations of 
H2 only and refer the reader to Naiita et al. ( 200?! ) for details on the inner disk 
gas component. 

As a symmetric molecule, H2 has no electric dipole moment and only magnetic 
quadrupole transitions are permitted, i.e. AJ = ±2. These lines lie in the mid- 
infrared, with energies characteristic of the temperatures at several AU in a disk, 
a potentially very interesting regime for understanding giant planet formation. 
Bitner et all (|2007l ) detected and spectrally resolved 3 lines (H2 S(l), S(2), S(4) 
corresponding to J = 3 — 1,4 — 2,6 — 4 respectively) in AB Aur. The linewidths 
locate the emitting region to 18 AU in radius but the rotational diagram indicates 
a temperature of 670 K, which is much hotter than the dust temperature at 
this radius. They infer that the gas is decoupled from the dust and excited 
by e ither protostellar UV or X-ra ys, or shocks from t he inf alling envelope (see 
also lBitner et al.ll2008l ). Similarly. iMartin-Zai'di et al.l (|2007l ) detect the H2 S(l) 
line in another intermediate mass Herbig Ae star, HD97048, and show that the 
high line-to-continuum ratio requires either an enhanced gas-to-dust ratio or an 
additional gas heating mechanism such as X-rays or shocks. 

In both cases, the inferred gas mass is very small. For AB Aur, where the 
multiple lines allow the level populations to be determined, = 0.5 ~ 

0.1% of the total disk mass. Effectively, only the top of the disk atmosphere is 
me asured due to the high optic al depth of the dus t at these wavelengths. Based on 
the lChiang fc GoldreichI (| 19971 ) two-layer model, ICarmona et al.l (j2008l ) calculate 
that mid-infrared H2 detections require gas temperatures twice as high as the 
dust and a gas-to-dust ratio greater than 1000. This is possible through dust 
settling to the mid plane but the low detection rates suggest th at such extreme 



2OIOI I. 



conditions are rare (jLahuis et al.l 120071 . iMartin-Zai'di et al 

Trace amounts, Mhi ^ 1 Mm , of atomic hydrogen should al so be present due to 
photodissociation of H2 ( Kamp. Freudling &: Chengalurl 200?! ) . Although beyond 
the sensitivity of current instruments, it will be possible to detect t he 21 cm line 
in many disks with the Square Kilometer Array ( Kamp et al.ll2008l ). 

Fortunately, there are many other currently observable species present in disks 
from ionized and atomic at the disk surface to molecular deeper within. An inter- 
esting case is the optical line of [O I] at 6300 A, which has b een detected in the sur- 

face o f strongly flared Herbig Ae/Be disks out to ~ 100 AU (jAcke. van den Ancker &: DuUemond 
20051 '). It originates from the photodissociation of OH which results in excited 



neutral atomic oxygen. 

Spitzer spectroscopy has provided many ne w results. 



_[Ne II 



II is one of the 

strongest lines and has an origin in YSO jets (jGiidel et al.ll2010l ). For sources 
without jets, the [Ne II] line is much weaker and traces hot gas on the disk 
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surface, possibly from a photoevaporate flow ( Pascucci &: Sterzik 20091 ). Ro- 
vibrational lines from the organic molecules, C2H2, HCN, and C O?, have also 



been d etected by IRS, in a bsorption through a n early-edge on disk bv lLahuis et al. 



and in emission by lCarr &: Najital (j2008l ). Assuming local thermodynamic 
equilibrium, both studies find molecular excitation temperatures in the range 
of 300 — TOOK, implying a location in the inner 3 — 6AU of the disks. Their 
production is due to a rich gas phase chemistry sustained by the sublimation of 
icy grain mantles in a warm molecular region abo ve the disk midpla ne, and is 

(Hooa) find a 



Pascucci et al 



strongly dependent on the stellar radiation field, 
deficit of HCN in disks around very low mass stars an d brown dwarfs, whi c h they 
attribute to a reduced photodissociation rate of N2. 'Pontoppidan eialJ (|201ol ) 
find almost no molecular emission toward 25 Herbig Ae/Be stars (the exception 
being one CO2 line). They suggest either photodissociation, lower gas-to-dust 
ratios, or masking by the stron g continuum as the ca use. 

The bulk of the hues in the ICarr k Najital (|2008l ) IRS spectrum of A A Tau 
are rotational lines of OH and H2O. These are important coolants and also of 
great astrobiological interest. Their derived high abundances require vertical 
and potenti a lly ra dial transport of ices from the inner, cooler regions of the disk. 



Salvk et al.l (120081) found similar re s ults t oward an additional two disks, and the 



large survey by Pontoppidan et al. ( 2010l ) shows that these lines are common in 
T Tauri disks. The complex origin of the gas-phase water and its implications 
for understanding dis k str ucture and evolution h ave been recently discussed by 
Meiierink et"all (j2009l ) and lCiesla fc Cuzzil ^200(i ). 



There are also many detectable molecular lines in the (sub-)millimeter wave- 
length regime. Observations of low energy rotational transitions, which can be 
excited even at the low temperatures and densities of the outer disk, provide the 
best constraints on the total gas mass in a disk. CO lines are the strongest on 
account of its high abundance and large dipole moment. The dominant source of 
opacity in this case is not the dust but the lines themselves. Interferometric CO 
images therefore show the radial gas temperature gradient, T oc r"*^, with values 
of q 0.4 — 0.7 for both T Tauri and H erbig Ae/Be stars (ICuilloteau fc Dutrey 



1998I . IPietu. Dutrev fc Guillot^ l2007l 'l. which is in agreement with the dust 



temperature derived from SED modeling (lAndrews &: WilhamsllioO^I ). Multi- 
transition observations trace different r = 1 s urface s providing information on 



the vertical temperature structure. iQi et al. show that the strong CO 



6-5 emission of the TW Hydra disk requires a hot surface layer with a gas-dust 
temperature difference of up to 30 K that can be explained by X-ray heating. Ob- 
servations of the rarer isotopologues, ^^CO and C^^O, probe deeper into the disk 
and show cooler midplane temperatures, ~ 13 K, and an order-of magnitude de- 



pletio n of CO relative to the interstellar medium (jPartois. Dutrev Guilloteau 
20031 ). This i s due to CO condensation ont o dust grains at a temperature 



of about 20 K (Sandford Sz Allamandola 



19931). Siniilarly high depletioii s have 



been found in other T Tauri disks (e.g., lOberg et al.ll20ld . lOi et al.l[20o3 ). The 
situation is quite different for disks around the more luminous Herbig AeBe 
stars where the dust is warmer than the CO condensation temperature through- 
out. Within the uncertainties of C O is otopologue a bund ance and dust opacity, 
Pietu. Guilloteau fc DutrevI (|2005l ) and lPanic et al.l (|2008l ) show that the gas-to- 



dust ratios in the AB Aur and HD 169142 disks are consistent with the interstellar 
medium value of 100. 

Rotational lines of other molecules further constrain the physical structure as 
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well as the chemistry of the outer disk. The first results were obtained using 
single-dish telescopes, measuring the globally averaged properties, and showed 
that other molecules including HCO"*" and H2CO were also strongly depleted as 
much as, if not more than, CO in classical T Tauri disks (ID utrev, Guillote au &: Guelin 
19971 : iThi. van Zadelhoff fc van Dishoeckl I2nn4l : Ivan Zadelhoff et al. 200 li). Pho- 



todissociated products such as CN were relatively abundant, however, which led 
to the picture of a photon dominated region at the surface of a warm intermedi- 
ate layer with T ~ 20 — 40 K, n(H2) ~ 10^ — 10^ cm~^, sitting on a cold, highly 
depleted midplane. The d etails in this picture are being filled in as new lines and 
s pecies are detected (e.g., Henning et al. I2OI0I ) . 

Thi. van Zadelhoff van Dishoeckl (|2004l l found a high deuterium fraction in 



TW Hydra, similar to values seen in dense cores an d attributed to cool tem- 
peratures and high CO depletion ( Caselli et al. 19991 ). Interferometric imaging 
resolves the D/H ratio, shows tha ,t it varies fro m ~ 0.01 to ~ 0.1 across the 
disk and peaks at ~ 70 AU radius ( Qi et al. 20081 ). a result with implications for 
understanding the high D/H ra tio in comets and the delivery of water to Earth 
( Ehrenfreund &: Charnlev 2000l ) . The suite of molecular line maps in Figure H] 
directly shows the decrease in DCO"*" 3-2 intensity at the center of the disk where 
the CO depletion is lower due to the higher temperatures. Deuterated species, in 
particular H2D"'", and other molecules with very low condensation temperatures 
such as N2H"'" should reveal the physical conditions and kinematics of the disk 
midplane in the outer disk. The lines are weak bu t H2D+ has been detected in 
one, possibly tw o, sources bvlCeccarelli et al.l (l2004l'l and N^H"*" has be en detected 
in several disks (|Dutrev et al.ll200i lOberg et al.ll20l"nl . lOi et al.ll2003l ). The cool 
temperatures of T Tauri star disks lead to C O freeze out and a richer chemistry 
that is absent in warmer Herbig AeBe disks (|6berg et al.ll20inl ). The reader is 
referred to BerginI ( 20091 ) for a more exhaustive discussion of disk chemistry. 

The full analysis of the many species and lines now detected in protoplanetary 
disks requires sophisticated models to handle the physical, th ermal, and chemical 
struc ture and the feedback of each of these on each other (e.g., Woitke. Kamp &: Thil 
2OO9I ). The complexity is such that, for T Tauri stars at least, there is no simple 
prescription for the gas mass. However, these observations provide substantial 
information about the disk structure and complement the broadband dust mea- 
surements. Disks have a dynamic atmosphere above a cool midplane: the tem- 
perature increases with height, the phase changes from molecular to atomic to 
ionized, and dust grains undergo sublimation and condensation of icy mantles as 
they move throughout. 

Finally, as this review is being written, the first results are coming in from 
Herschel and its far-infrared spectrometers. The [ 01] 63 line is gener ally the 
strongest line in this regime and widely detected ( Mathews et al.l 2010l ). Com- 
bined with gr ound based CO measurements, it provides strong constraints on 
the gas mass ( Thi et al. 2010l ). Multiple transitions of this and other species 
can be detected in some sources and allow the excitation conditions of the gas 
to be determined ( Sturm et al.l 2010l ). However, the weakness of H2O lines at 
sub-millimeter wavelengths indicates that almost all the water is fr ozen onto icy 



mant les around large dust grains that have settled to the midplane ([Bergin et al. 
20ld ). In HerscheVs 3 year lifetime, we expect to learn much more about the 



gaseous content of disks in the giant planet-forming zone. 
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4.5 Dependence on stellar mass 

Disks are observed around a wide range of stars from very low mass to inter- 
mediate mass Herbig Ae/Be stars. Infrared excesses from disks around b rown 
dwarfs are also detected at a similar frequency to stars ( Luhman et al. 20051 ) and 
masses for a handful have been measured fr om millimeter wavelength measure- 
ments pcholz. Jayawardhana &: Woodll200()l ). The relationship between disk and 
stellar properties informs models of disk formation and evolution. 

Higher mass stars require more material to pass through a disk and we might 
there fore expect to see a positive corre lation between disk mass and stellar mass 

(e.g., Natta. Grinin &: Manniiigi 2000l ) . Figure [SI c ompiles (sub-)mi l limeter mea- 

surements of Class II YSO disks from the literature ( Acke et al. I2OO4 lAlonso-Albi et al 



20091 : [Mannings fc Sargentl[2ciool : iNatta. Grinin fc Manningsll2000l : IScholz. Javawardhana fc Wood 



20061 ) and shows that there is a large scatter, ~ 0.5 dex, but confirms that disk 
masses tend to be lower around low mass stars such that the ratio, Ma/M^ ~ 0.01. 

The relationship breaks down for the most massive stars. There is no (sub- 
)millimeter detection of a disk around an optically visible O star. The limits 
on the sensitive Orion proplyd survey provide the most stringent constraints to 
date,, Md/M^ < 10"^ for > WMq (|Mann fc Wilhamsl l2009bl v This may be 
due to very high photoevaporation rates so that disks are not detectable by the 
time an O st ar is optically visible or d ue to an altogether different star formation 
mechanism ( Zinnecker &: Yorkell2007l l. 

Disk size measureme nts, whether from optical observations of O rion proplyds 
(IVicent e &: Alvesll2005l ) or millimeter visibilities of Ophiuchus disks ( Andrews et al. 
I2OO9. . .2010131 ) , are of a restricted range of stellar masses, M^, 0.3 — 2 Mq , and 



do not follow a clear trend. Moreover, iPanic &: Hog erheiidd (|2009l ) argue that 
observations of Herbig Ae/Be stars at millimeter wavelengths have been biased 
toward large, bright disks. ALMA will soon provide the resolution and sensitivity 
required to measure disk radii in a more representative sample. 

Disk stru cture does appear to be dependent on stellar mass, at least for low 
mass stars. Szilcs et al. ( 20ld ) find that relatively blue mid-infrared colors imply 
smaller scale heights for disks around very low mass stars than low mass stars 
in the same star-forming region. It is harder to make a direct comparison to 
higher mass Herbig Ae/Be stars in other regions because of the difficulty in 
distinguishing initial conditions from evolution due to grain growth and settling. 



4.6 Section summary 

• Disk masses are best measured from (sub-)millimeter observations of ther- 
mal continuum emission from dust. Substantial, and uncertain, corrections 
are required for the grain size distribution and gas-to-dust ratio. 

• Assuming an ISM gas-to-dust ratio of 100 and ignoring the mass in bodies 
larger than 1 mm, the median mass of protoplanetary disks around Class II 
YSOs with spectral types GKM is 5Mjup. 

• Protoplanetary disks have power law surface density profiles with an expo- 
nential taper that produces a soft edge anywhere from 20 to 200 AU from 
the star. The power law index is approximately —1 but is extrapolated and 
highly uncertain within the central ~ 20 AU. 

• Surface densities are generally below the threshold for gravitational insta- 
bility and velocity profiles are Keplerian. 
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• Mid-infrared through milhmeter spectroscopy reveals a warm molecular 
layer with PDR-like chemistry that is strongly dependent on the stellar 
luminosity and a cold midplane depleted of molecules. 

• From brown dwarfs to B stars, disk masses scale with the stellar mass. The 
median ratio is about 1% but there is a large dispersion, ±0.5 dex. 



5 DISK LIFETIMES 



One of the most fundamental parameters on disk evolution studies is the lifetime 
of the disk itself. This is not only because it reflects the relevant time scale of 
the physical processes driving the dissipation of the disk, but also because it sets 
a limit on the time available for planet formation. 

Even though the mass of the gas initially dominates that of the dust by two 
orders of magnitude in primordial disks, the dust is much easier to observe and 
thus most of the constraints on the lifetimes of circumstellar disks have been 
obtained by observing the thermal emission of the dust grains. These particles 
absorb stellar light and reradiate mostly in the 1 fim to 1 mm range. Because 
the temperature of the dust decreases with the distance from the central star, 
different wavelengths trace different disk radii (for a given stellar luminosity). 
In what follows, we summarize the observational constraints provided by near- 
infrared and mid-infrared observations, mostly from Spitzer, which has performed 
the most sensitive and comprehensive surveys of disks in star-forming regions to 
date. Toward the end of the section, we also briefly discuss current constraints 
on gas dispersal time scales. 



5.1 Near-infrared results: the inner disk 



Because there is a very well established correlation between the presence of near- 
infrared excess (1 — 5 um) and the occurrence of spectroscopic signatures of ac- 
cretion ( Hartigan. Edwards &: Ghandourlll995l ). it is possible to investigate the 
lifetime of inner accretion disks (R ^ 0.1 AU) by studying the fraction of stars 
with near-infrared excess as a function of stellar age. Early studies of nearby 
star-forming regions found that 60 — 80% of stars younger than 1 Myr present 
measurable near- in frared excesses, and that no more than 10% of stars older than 
10 Myr do so (e.g., IStrom et al] 1 19891 ). Given the large uncertainties associated 
with model-derived stellar ages, it has been argued that individual star- forming 
regions lack the intrinsic age spread necessary to investigate disk lifetimes from 
individually derived ages an d the apparent a ge dispersion is mostly driven by the 
observational uncertainties (lHartmannll200lh . However, sir nilar studies based on 
the disk frequencies in c lusters with a range of mean ages ( Haisch. Lada Sz Ladal 
200 ll : iHillenbrandl I2OO5I ) have led to essentially the same result: the frequency of 
in ner accretion disks steadily decreases from < 1 to ~ 10 Myr (e.g., see Figure 2 
in lWvattJ (I2OO8I II. 

The inner disk fractions observed in both young stellar clusters and in the 
distributed population of pre-main sequence stars in star-forming regions are 
consistent with median disk lifetimes of between 2 and 3 Myr. The decrease 



has a lso been modeled as an exponential with e-folding time 2.5 Myr (iMamajek 
2OO9I ). The uncertainty in this number is mainly due to the difficulty in deter- 



mining stellar ages accurately at these early times, and the uncertain duration of 
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star formation in a cluster. There is also a wide dispersion in the lifetimes: some 
objects lose their inner disks at a very early age ( ^ 1 Myr), even before they be- 
come optically revealed and can be accurately placed in the Hertzsprung-Russell 
diagram, whereas other objects retain their accretion disks for up to 10 Myr. 



5.2 Spitzer results: the planet-forming regions of the disk 

Because only circumstellar dust very close to the star (i?^0.1AU) becomes 
hot enough to be detectable in the near- infrared, observations at these wave- 
lengths provide no information on the presence of circumstellar material beyond 
~ 0.5 AU. As a result, disk lifetime studies based on near-infrared excesses always 
left room for the possibility that stars without near-infrared excess had longer- 
lived outer disks with enough material to form planets at radii that could not be 
detected at these short wavelengths. 

IRAS and ISO both had the appropriate wavelength range to study the planet- 
forming regions of the disk (i? ~ 0.5 — 20 AU) but lacked the sensitivity needed 
to detect all but the strongest mid- and far-infrared excesses in low-mass stars 
at the distances of the nearest star-forming regions. During its 6-year cryogenic 
mission, Spitzer provided, for the first time, the wavelength coverage and the 
sensitivity needed to detect very small amounts of dust in the planet-forming 
regions of thousands of YSOs. The wealth of data made available by Spitzer has 
not only firmly established the dissipation time scale of primordial disks, but also 
made it possible to address second-order questions such as the effect of stellar 
mass, multiplicity, and external environment on disk lifetimes. 
5.2.1 Stellar CLUSTERS AND ASSOCIATIONS Just as was done for the inner 
disk with near- infrared studies, the dissipation time scale of regions farther out 
in the disk can be investigated by observing the fraction of stars that show mid- 
infrared excesses in clusters of different ages. Reaching the stellar photospheres 
of all the targets is needed to unambiguously establish the fraction of them that 
have an infrared excess indicating the presence of a disk. However, while at 
8.0 //m Spitzer was able to detect the stellar photospheres of solar-type stars 
at distances of up to 1 kpc, at longer wavelengths it was only sensitive enough 
to detect solar-type photospheres within ~ 200 pc. Because there are very few 
stellar clusters within 200 pc, we focus on the results from IRAC (3.6 to 8.0 /^m) 
observations. 

Comparing the disk fractions of stellar clusters presented by different Spitzer 
studies is not completely straightforward as the methods used to decide which 
objects are members of the cluster, the disk identification criteria, sensitivities, 
and ranges of stellar masses considered vary from study to study. 

Very young embedded clusters (age < 1 Myr) such as Serp ens an d NGC 1333 



show disk fractions of the order of 70 to 80% ( Gutermuth et al . 2008, Winston et al 



20071 ). This implies that pre-main sequence stars without a disk, showing pho- 



tospheric IRAC fiuxes, are seen even in extremely young star-forming regions 
Clusters in the 2 - 3 Myr range, such as IC 348 a.nd NGC 2264, show IRAC disk 
fract ions of the order of 40 to 50% ( Lada et apbood . ISung. Stauffer fc Bessell 



20091 ). whereas regions with estimated ages around ~ 5 Myr, such as Upper 



Scorpius and NGC 2362 exhibit disk fractions that are already below 20%. By 
~ 8 — 10 Myr old, primordial disks with IRAC excesses become exceedingly rare 
as attested by their very low inciden ce ( ^ 5%) in regions such as the TW Hydra, 



£7 Ori, and NGC 7160 associations (iHernandez et al.l 120071 . ISicilia-Aguilar et al 
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The main results on disk lifetimes from IRAC studies of young stellar clus- 
ters listed above are indistinguishable from those obtained by pie-Spitzer near- 
infrared studies. This is not too surprising considering that IRAC wavelengths 
only trace regions inward of ~ 5 AU around solar type stars. These regions are 
further out in the disk than those traced by near-infrared observations, but still 
exclude the bulk of the disk material. Thus, the question remains: can primor- 
dial outer disks {R > 5AU) survive beyond lOMyr? To address this question, 
we must turn to the longer wavelength camera, MIPS, on Spitzer, and studies of 
star-forming regions within 200 pc, where 24 ^um observations are sensitive enough 
to reach the stellar photosphere of their targets and unambiguously identify the 
presence of a disk. 

5.2.2 The DISTRIBUTED POPULATION OF PRE-MAIN SEQUENCE STARS Firmly 

establishing the longevity of primordial disks was one of the central goals of two 
Spitzer Legacy P rograms: "From Molecular Cores to Planet-forming Disks" (cd2: 



Evans et al.ll2003i) and "Formation and Evolution of Planetary Systems" (FEPS: 



Mever et al.ll200()l ). 



As part of the c2d project, Spitzer observed over 150 WTTS associated with 
the Chameleon, Lupus, Ophiuchus, and Taurus star-forming regions, all within 
200 pc of the sun. As one of the main goals was to establish whether the outer disk 
could significantly outlive the inner disk, the c2d sample of WTTS was mostly 
composed of X-ray identified and spectroscopically confirmed pre-main sequence 

stars without evidence for an i nner accretion dis k. 

Results from the c2d project (ICieza et al ] l20n7l . IPadgett et J]l200fil . IWahhai et ^. 



show that ~ 80% of young WTTS present 24 nm fluxes consistent with 



bare stellar photospheres. Because 24 /im observations are sensitive to very small 
amounts of micron sized dust (<C 1 M®) out to tens of AU from the central star, 
this suggests that once accretion stops and the inner disk dissipates, the entire 
disk goes away very quickly. The fraction of c2d WTTS with 24 /xm excesses 
is a function of stellar age. Approximately 50% of the WTTS younger than 
^1 — 2 Myr do not have a disk, suggesting that a significant fraction of the 
entire pre-main sequence population lose their disks by an age of ~ 1 Myr. How- 
ever, none of the WTTS in the c2d sample older than ~ 10 Myr has a detectable 
disk. Also, many of the WTTS disks have very low fractional disk luminosities, 
(-^disk/-^* < 10~^), and are thus more consistent with optically thin debris disks 
than with primordial disks. 

Very similar results were obtained by the FEPS project, which observed 314 
solar-type stars with a median distance of 50 pc and ages ranging from 3 Myr 
to 3 Gyr. The youngest age bin in the FEPS study included 34 targets younger 
than 10 Myr and a mean age of ~ 5 Myr. They find that only 4 out of those 
34 targets had optically thick primordial disks, whereas 5 of them had optically 
thin debris disks. The second youngest bin in their study had 49 targets with 
estimated ages in the 10 to 30 Myr ra nge. In this age bin, they identified 9 
debris disk and only one primordial disk ( Carpenter et al. 20091 ). The only target 



with a primordial disk in this age range was PDS 66, a member of the Lower 
Centaurus Crux association with an estimated, but uncertain, age of 12 Myr 
(jPreibisch fc MamaieklboOSh . 



Both the c2d and FEPS results strongly suggest that, while there is a wide 
dispersion on disk lifetimes, ~ 10 Myr is a firm upper limit for the longevity 
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of primordial circumstellar disks around solar-type stars. They also show that 
there is a significant overlap in the age distributions of primordial and debris 
disks. Exactly how and when protoplanetary disks evolve into planetary debris 
disks remains an open question and is discussed in §6.4[ 



5.3 Dissipation timescale 

The fact that very few objects lacking near-infrared excess show mid- infrared ex- 
cess emission implies that, once accretion stops and the inner disk clears out, the 
entire disk dissipates very rapidly. Th i s is also supported by (sub-) millime ter ob- 
servations ( Andre Sz Montmerlj 1994 . Andrews Williams 2005 . 2007al ) show- 
ing a very strong correlation between the detectability of a disk at these long 
wavelengths and presence of an inner accretion disk. This implies that the vast 
majority of pre-main sequence stars in any given population are either accreting 
CTTS with excess emission extending all the way from the near-infrared to the 
sub-millimeter or have bare stellar photospheres. Based on the small incidence 
of objects lacking an inner disk that have evidence of an outer disk, the dissipa- 
tion timescale of the entire primordial disk once accretion stops has repeatedl 



been estimated to b e ^O.SMyr (e.g.. ICieza et al.M2007l . ISkrutskie et al.Ml99i 



I 



Wolk fc Waited Il996l ). This finding, that circumstellar material can survive at 



all radii {R ^ 0.1 — 200 AU) for several Myr, while later the dissipation of the 
entire disk occurs in a much shorter timescale is known as the "two-time-scale" 
problem. D isk evolution models combining vis cous accretion with photoevapo- 
ration (e.g., Alexander. Clarke fc Pringle 2006al lbl) have successfully been able to 
reproduce this behavior (see ^6.2p . 



5.4 Dependence on stellar mass 

Young stellar clusters and associations provide an excellent opportunity to inves- 
tigate disk lifetimes as a function of stellar mass as they contain large samples of 
mostly coeval stars spanning a large range of stellar masses. In the ir study of the 
5 Myr old Upper Scorpius OB association. Carpenter et al. ( 20061 ) obtained 4.5, 
8.0 and 16 //m photometry of 204 members with masses ranging from ~ 0.1 to 
20 Mq. They find that while ~ 20% of their 127 K and M-type targets (masses 
~ 0.1 — 1.2 M0) are surrounded by optically thick, primordial circumstellar disks, 
none of their 30 F and G-type stars (masses ~ 1.2 — 1.8 Mq) had any evidence 
for a disk at wavelengths < 16 fj,m. 

An almost identical result was obtained by Dahm Sz Hillenbrand ( 200?! ) in the 
also 5 Myr old cluster NGC 2362. They found a total IRAC disk fraction of 
~ 20% (7% of strong excesses and 12% of weaker excesses) among 220 stars 
with estimated masses below 1.2 Mq and a complete absence of IRAC excesses 
among 33 s tars with estimate d masses above 1.2 Mq. MIPS observations of 
NGC 2362 (ICurrie et al.l l2009l) and other 5 Myr old clusters such as A Orionis 
(iHernandez et al.l l2009l )" do show 24 fxm excesses for objects with masses above 
1.2 Mq, but they are all consistent with the optically thin emission expected 
from debris disks. Overall, the Spitzer results suggest that, though primordial 
circumstellar disks can last for up to 10 Myr around solar and lower-mass stars, 
disk lifetimes are a factor of ~ 2 shorter around higher mass objects. The higher 
accretion rates and radiation environmei it are likely to be responsible for the 
shorter disk lifetimes in higher mass stars (jCalvet et al.ll2005l . lGarcia Lopez et al. 
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20061 . lHillenbrandl2008l l. 



At the lower end of the mass distribution, the statistics are much poorer. The 
disk fractions of brown dwarfs in relatively young regions (age 1 — 3 Myr) such as 
Taurus, Chamaeleon I, and IC 348 has been fou nd to be between 40 and 50% with 
statistical uncertainties as large as 10 to 20% fiGuieu et al. 2007 . iLuhman et al. 



200,1 iMonin et all l20inl ). For the 5 Myr Up per Scorpius region, the reported 
values range from 37 ± 9% (jScholz et al.ll2007l ^ to lll^ 3% (iRiaz. Lodieu fc Gizia 
2OO9I ). Also, 3 of the 5 known brown dwarfs in the 10 Myr old T W Hydrae 



assoc iation have infrared excesses indicating the presence of a disk (|Riaz &: Gizis 
20081 ). Taken together, these observations suggest that the dissipation time scales 



of disks around substellar objects are at least as long as those of solar-mass stars, 
although significantly longer disk lifetimes cannot be ruled out based on the 
available data. 



5.5 Gas dispersal 



CO observations toward disks that are warm enough to pre vent freeze-out hint at 
a decrease in the gas-to-dust ratio with evolutionary class (jChapillon et al.ll20ld : 
Dutrev. Guilloteau &: Simon l2003l ^. Unlike dust, however, there is no generic 
tracer of the gas mass in young disks (see ^4.4.2p and systematic studies of gas 
dispersal with stellar age has been studied by observing accretion indicators. 

Using t he equivalent widths and velocity profiles of the Ha line to identify 
accretors, Fedele et al. ( 2010l ) recently studied the fractions of accreting objects 
in stellar clusters with ages in the 1—50 Myr range. They find no accreting objects 
in the clusters older than 10 Myr down to an accretion sensitivity estimated to be 
^ 10""^^ M0yr~^. They also find that, in most clusters, the fraction of accreting 
stars is systematically lower than the fraction of objects with Spitzer excesses 
in IRAC bands. This is expected as ~ 20% of non-accre t ing pre-main sequence 
stars (i.e., WTTS) have IRAC excesses (jCieza et al.l 120071 . iDamjanov et al.1 120071 . 
Wahhaj et aDl20inl ). while the fraction of accreting pre-main sequence stars (i.e., 
CTT S) lacking IRAC excesses is significantly smaller, of the order of ^ 2-5% 
(e.g., ICieza et al.ll2010l . iMuzerolle etnill20inl i. 

Accretion only shows the presence of gas in the inner disk. Consequently, our 
current understanding of the time scale of the gas dispersal is in an analogous 
situation to that of the dust dissipation prior to Spitzer: although the longevity 
of gaseous inner disks is well established to be ^ 10 Myr, the possibility of sig- 
nificant amounts of gas remaining for longer period of times at larger radii still 
exists. However, according to the "UV-switch" model (Alexander et al., 2006), 
photoevaporation removes all circumstellar gas very quickly (<C 1 Myr) once ac- 
cretion stops (see §6.2p . Herschel studies currently underway of sensitive gas 
tracers, such as the 63.2 /im [O I] line, should be able to test this prediction. 



5.6 Environmental influences 

Protoplanetary disks are detected in a range of environments from low mass, 
sparsely populated molecular clouds to massive, dense stellar clusters. In nearby, 
well studied, low mass star-forming regions, there is a remarkable similarity in the 
averag e disk properties from region to region as demonst rated from mid-infrared 
colors ( Fang et al. 20091 ) and spectra ( Furlan et al. 2009) and also in thei r mass 
distributions and sub- millimeter colors ( Andrews &: Williams! 2005 , 2007al ) . This 
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is not unexpected given that disk sizes are orders of magnitude smaller than 
typical star-to-star distances except perhaps in the densest parts of very young 
protoclusters and disrupti ve encounters with a passing star is not a common 
evolutionary determinant (jScally &: Clarkell2nnih . Gravitational perturbations 
have a strong impact on disk evolution, however, in close binary or multiple 
systems. 

The median disk lifetime of 2-3 Myr (discussed in ^ is derived from surveys of 
many clusters. There is a rather small dispersion in the disk fraction at any given 
age despite the different cluster sizes and whether or not they contain massive 
stars. These mid-infrared observations are unable to measure the properties of 
the outer parts of the disk, however, which are more susceptible to external 
influences. Here, we discuss the influences of binary stars, massive stars, and also 
metallicity on disk lifetimes. 

5.6.1 Dynamical Disruption in Binari es About half of fie ld stars are 
in binary or higher order multiple systems (Ragh avan et al. 2010l ). The or- 
bital resonances in th e se sys tems have a profound influence on disk evolution. 
Artymowicz &: Lubow ( 19941 ) show that (coplanar) disks around each star in a 



binary system are truncated at the outer edge, and a circumbinary ring about 
both stars are truncated at the inner edge. As a rough guide for near circular 
orbits with semi-major axis a, the circumprimary disk is limited in size to ~ a/2 
and the circumbinary disk's inner edge is ~ 2a. Higher eccentricities lead to 
greater disk erosion. 

There are a few resolved disk images that confirm this general pattern. GG Tau 
is perha ps one of the best studied circurn binary disks with an inner radius of 
180 AU ( Guilloteau. Dutrey &: Simo ^ Il999h although it a ppears to be bigger tha n 
expected given the a = _32AU,e = 0.34 binar v orbit (iBeu st fc DutrevI lioOfil'). 



gr 

UY Aur (jDuvert et al.lll998l ) and CoKu Tau 4 (jlreland fc Kr aus 2008) are ad- 



ditional examples of primordial circumbinary disks but they appear to be rare 
exceptions rather than the rule. Outwardly truncated circumstellar disks, only a 
few te ns of AU in diameter , have also been detected in binary systems, L1551- 
IRS5 (iRodriguez et al.l [l998l'l . GQ Lup (iDai et al.ll2o"Tol ^. HD 98800 and Hen 3-600 
(|Andrews et al.ll2010al ). 

High resolution speckle imaging studies searching for a connection between 
binaries and premature disk dissipation initia lly yielded mixed and inconclusive 
results. Ghez. Neugebauer &: Matthews ( 19931 ) surveyed ~ 70 pre-main sequence 
stars in Taurus and Ophiuchus and concluded that the incidence of close binaries 
(o < 50 AU) in WTTS is enhanced with respect to that of C TTS. This was not 
confi r med, however , by s ubsequent larger studies in Taurus (|Kohler &: Leinert 
19981 . iLeinert et aDll993ll. Yet a survey of ~ 160 pre-main sequence stars in 



Ophiuchus by Ratzka. Kohler fc Leinert ( 20051 ) showed that YSOs with infrared 
excesses tend to have fewer companions and are at smaller projected separations 
than diskless YSOs. 

Here, again, the uniformity and sensitivity of Spit zer studies have improved 
our understanding of the situation. ICieza et al. (|2009l ) combined the results from 
several multiplicity surveys of pre-main sequence stars with Spitzer data of four 
star-forming regions and showed that the distribution of projected separations of 
systems with mid-infrared excesses are in fact significantly different from that of 
systems without. Binaries with projected separations less than 40 AU are half as 
likely to possess a disk than those with projected separations in the 40 — 400 AU 
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range. iDuchend (|201(3l l finds a similar result with a somewhat smaller sample. 



(Sub-)millimeter fluxes are also known to be lower for similarly close binaries 
indicating that the bulk of the disk has been lost (lAndrews k Wilhain] l2005l. 
I.Tensen. Mathieu fc Full^ll994l . Il99fil . lOsterloh fc Beckwithlll9'95f r 

Even though several factors (e.g., the incompleteness of the census of close bina- 
ries, the use of unresolved disk indicators, and projection effects) tend to weaken 
its observable signature, the effect of multiplicity on disk lifetimes could be very 

strong. The distribution of physical separations, a, in solar- ty pe pre-main se- 

queiice binaries is expected to peak around 30 AU as in field stars (jDuquennov &: Mayor 
199ll . lRaghavan et al.l[201ol l. Hence the disks around individual stars in most bi- 
nary systems should have tru ncation radii of the order of (0.3 — 0.5)a ~ 10 — 15 AU 
( Papaloizou &: Pringld T977I ) . These truncation radii are about an order of mag- 
nitude smaller than the typical radii of disks around single stars (see ^4.2p . For 
an accretion disk with 7 = 1 (see ^4.3.ip . the viscous time scales linearly with 
radius and is therefore also about an order of magnitude smaller. This implies 
that the lifetimes of disks around the individual components of most binary sys- 
tems should be ~ 10% of those of single stars or about 0.3 Myr. These very short 
lifetimes for the stars in medium-separation binary systems may explain both the 
presence of very young diskless stars and the large dispersion in disk dissipation 
time scales. 

5.6.2 Photoevaporation by Massive S tars Most stars form in large 
clusters with hundreds if not thousands of stars ( Lada &: Ladall2003l . lMcKee fc Wilhamsl 
1993). Such large stellar groups are highly likely to contain an O star which bathes 



neighboring stars and their disks in a UV radiation field that may be several or- 
ders of magnitude above the interstellar average. The effect is to rapidly erode 
the more loosely bound outer parts of a disk. 

The spectacul ar Hubble images of photoevaporatin g disks embedded in small 
ionized cocoons (|Ballv. O'Dell fc McCaughreanll200nl ) make a stunning contrast 



to th e large Taurus disks in their more benign low radiation environment (jPadgett et al, 
I999I ). The pictures are somewhat misleading, however, in that the central re- 



gions of Orion disks are not strongly affected. Except for disks very close to the 
O stars where ionizing photons directly impinge on the disk surface, the disk 
develops a thick photon dominated region with a temperature ~ 10^ K. The cor- 
responding thermal velo city, ^ 3 km s~^, im plies that gas is bound to the central 
star for radii ^ 100 AU (j Adams et al.ll2004l ). Gas pressure can cause mass loss at 
smaller radii b ut the evapor ation timescale within solar system scales, < 50 AU, 
is tens of Myr (IClarkell2007l ). 

The comparison of disk masses in Orion with Taurus and Ophiuchus in Fig- 
ure [1] illustrate these points. There is a deficit of massive disks at the center 
of the Trapezium Cluster but such disks would likely have been the largest and 
therefore the most susceptible to photoevaporation by the 06 star, 6^ OriC. The 
discrepancy at the upper end of the mass distribution is roughly consistent with 
the smaller median size of Orion proplyds compared to Taurus disks. Further, 
the fraction of disks with at least a MMSN, lOMjup, within 60 AU is similar, 
~ 11 — 13%, in Taurus, Ophiuchus, and Orion and, i ncidentally, comparable t o 
the detection rate of Jupiter mass extrasolar planets (jMann &: Winiamsl[2009bl ;). 

More massive disks are found further away from the center of the Trapezium 
Cluster (lEisner et al.ll2008l . iMann fc Winiamsll2009al l and the disk mass distribu- 
tion of the full Orion region more closely resembles those in Taurus and Ophiuchus 
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(see Figure [T|) . The large difference in mass between the most massive disks in 
and out of the cluster center comp ared to the difference in stellar de nsity rules 
out stellar encounters as the cause ( Olczak. Pfalzner &: Spurzem 20061 ). 

VLA and SMA observations do not currently have sufficient sensitivity to study 
more distant regions. However, Hubble and Spitzer observations find morpho- 
logical features_similai^ massive star-forin i ng re - 
gions (iBalog et al. l l200fil : ISmith. Ballv fc Morsel l2003l ;i . iHernandez et al.l (|2008l ) 
find that the 7 Velorum cluster has a low fraction of sources with infrared ex- 
cesses for its age and consider disk photoevaporation as a possibility but also note 
the large uncertainty in the cluster age. The disk fractio n decreases by ab out a 
factor of 2 within the central 0.5 pc of the Rosette nebula tealog et al. 200?! ) and 
the central 1 pc of S Monoceros (jSung. Stauffer &: BesselL,2009. ) but the statistics 
are limited by small numbers. External evaporation can only remove the inner 
disk on such short timescales with v ery strong radiation fields that di rectly ionize 
the disk and boil it off at ~ 10^ K (|Adams et al.ll2004l . IClarkell2007l l. It is hard 
to reconcile this with a parsec-scale s phere of influence of a mass ive star unless 
t he stellar orb i ts are highly eccentric ( Storzer &: Hollenbach 19991 ). 

Balog et al. ( 20081 ) show that the Spitzer 24 /im dusty tails in these regions are 



not detected in Paa and are effectively gas-free. They suggest that photoevapo- 
ration removes the gas very quickly but that a dus ty reservoir is replenish ed by 
the collisions of large bodies left behind in the disk. iThroop fc Bdlvl (|2005l l con- 
sider the rapid enhancement of the dust-to-gas ratio in photoevaporating disks 
as a potential trigger for planetesimal formation. External evaporation has also 
been considered as a possible explanati on of the steep drop in the su rface density 
of Kuiper Belt Objects beyond 50 AU (I.Tewitt. Luu fc Trujillolf]^ ). 

5.6.3 The effect of metallicity The effect of the metalhcity, or the 
initial dust-to- gas content, on the evolution of protoplanetary disks has not 
been well studied, at least at the early stages. It is known that the metal- 
l icity of the host star i s strongly correlated with the presence of hot Jupiters 
(IValenti fc F ischer^ l 200,5l). This may be due to a n enhanced rate of planetesimal 
formation (jjohanse n. Youdin &: Mac Low 2009) but lon ger disk lifetimes may 
also play a role. In a small sample, D'Orazi et al. ( 20091 ) do not find any signif- 
icant differencej^n the metallicities of Taurus and Orion WTTS or CTTS. How- 
Yasui et al. (|2009l l showed the the inner disk fraction, as measured by near- 



ever, 



infrared exc esses, is sign i ficant ly smaller in low metallicity clusters at the edge of 
the Galaxy. Yasui et al. ( 20ld ) extend this work to other distant, low metallicity 
clusters of different ages and find a median disk lifetime, < 1 Myr. They suggest a 
strong dependence of the lifetime on metallicity, which is in qualitative agr eement 
with photoevaporative models of disk evolution ( Ercolano &: Clarke! 20ld ). 



5.7 Section summary 

• Disk lifetimes have a median between 2 and 3 Myr but vary from less than 
1 Myr to a maximum of 10 Myr. Multi-wavelength observations show that 
the dissipation timescale is very similar at all radii. 

• Circumstellar disks dissipate faster around high-mass stars than around 
solar mass stars. Very low mass-stars and brown dwarfs have disk lifetimes 
at least as long as, and likely longer than, those of solar mass stars. 

• Disk lifetimes are significantly shorter around medium separation binaries 
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(a ~ 5 - 100 AU). 

• Photoevaporation from massive stars erodes the outer parts of disks but 
generally leaves their interiors ( ^ 50 AU) intact. The Hfetimes and masses 
are only significantly reduced for disks that lie within a few tenths of a 
parsec from an O star. 

• Disks in low metallicity environments at the edge of the Galaxy have me- 
dian lifetimes < 1 Myr, which are substantially shorter than in the Solar 
neighborhood. 



6 DISK EVOLUTION 



Understanding the physical processes that drive the evolution of primordial cir- 
cumstellar disks, as they evolve from optically thick to optically thin, is crucial for 
our understanding of planet formation. Disks evolve through various processes, 
including viscous accretion, dust settling and coagulation, dynamical interac- 
tions with (sub-)stellar companions and forming planets, and photo-evaporation 
by ultraviolet and X-ray radiation. In this section we summarize the models and 
observational constraints for the different processes that control the evolution of 
primordial circumstellar disks. 



6.1 Viscous transport 

To first order, the evolution of primordial disks is driven by viscous transport. 
The accretion from the inner disk onto the star is relatively well understood 
and well constrained observationally. The large velocity widths, intensities, and 
profiles of emission lines such as Ha, Br7 and Ca II can be successfully repro duced 
by magnetospheric accretion models ( Muzerolle. Calvet &: Hartmann 19981 ). The 
ph ysical mec h anism s that drive the radial transport across the disk are reviewed 
bv lArmitagj toid ). 

Viscous evolution models are broadly consistent with the observational con- 
straints for disk masses and sizes, and the de crease in accretion rate over time 
(iHartmann et al.lll998l . bueso fc GuillotJ 1200,^1 ^. ^lowever, pure viscous evolution 
models also predict a smooth, power-law evolution of the disk properties. This 
secular disk evolution is inconsistent with the very rapid disk dissipation that 
usually occurs after a much longer disk lifetime (i.e. the "two-time-scale" prob- 
lem). Viscous evolution models also fail to explain the variety of SEDs observed 
in the transition objects discussed in ^ These important limitations of the vis- 
cous evolution models show that they are in fact just a first-order approximation 
of a much more complex evolution involving several other important physical 
processes. 



6.2 Photoevaporation by radiation from the central star 

Together with viscous accretion, photoevaporation is one of the main mechanisms 
through which primordial circumstellar disks are believed to lose mass and even- 
tually dissipate. Photoevaporation can be driven by energetic photons in the 
far-ultraviolet (FUV: 6eY < hi^ < 13.6 eV), extreme-ultraviolet (EUV: 13.6 eV 
< hv < 0.1 keV) and X-ray {hv > 0.1 keV) energy range. Photons in each en- 
ergy range affect the disks in different ways, and the relative importance of FUV, 
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EUV, and X-ray photoevaporation is still not well understood. Photoevaporat- 
ing photons can originate both from nearby massive stars and from the central 
star itself. The former scenario was discussed in ^5.6.21 and here we focus on the 
latter. 

Ea rly photoevaporation ni odels focused on the effect of ionizing EUV radiation 
(e.g., iHollenbach et al.lll994l ) on circumstellar gas around early-type stars. Later 
disk evolution models, known as "UV-swit ch" models, combine viscous evoluti on 
with photoevaporation by EUV ph otons (jAlexander. Clarke &: Pring i3 l2nn6al lbl: 



Clarke. Gendrin &: Sotomavorll200ll ) to tackle the "two-time-scale" problem of T 



Tauri evolution (i.e., the sudden dispersion of the entire disk after much longer 
disk lifetimes, see ^5.3p . According to these models, extreme ultraviolet (EUV) 
photons originating at the stellar chromospheres of low-mass stars ionize and 
heat the circumstellar hydrogen to ~ 10^ K. Beyond a critical radius, ~ 10 AU 
for solar mass stars, the thermal velocity of the ionized hydrogen exceeds its 
escape velocity and the material is lost in the form of a wind. 

At early stages in the evolution of the disk, the accretion rate dominates over 
the evaporation rate and the disk undergoes standard viscous evolution: material 
from the inner disk is accreted onto the star, while the outer disk behaves as a 
reservoir that resupplies the inner disk, spreading as angular momentum is trans- 
ported outwards. Later on, as the accretion rate drops to the photoevaporation 
rate, ~ 10"^'' — 10~^ Mq yr~^ in the models, the outer disk is no longer able to 
resupply the inner disk with material. At this point, the inner disk drains on a 
viscous timescale ( ^ 10^ yr) and an inner hole of a few AU in radius is formed 
in the disk. The inner disk edge is now directly exposed to the EUV radiation 
and the disk rapidly photoevaporates from the inside out. Thus, the UV-switch 
model naturally accounts for the lifetimes and dissipation timescales of disks as 
well as for SEDs of some pre-main sequence stars suggesting the presence of large 
inner holes. 

Recent photoe vaporation models include X -ray ('Owen et al. 2O10l. e.g..) and/or 
FUV irradiation (ICorti. Dullemond fc Hollenb ach 2009: Gorti fc Hollenbachll2nn9l . 
e.g.,) in addition to the EUV photons. Because X-rays and FUV photons are 
able to penetrate much larger columns of neutral gas than EUV photons, they 
are able to heat gas that is located both deeper in the disk and at larger radii. 
Thus, while EUV induced photoevaporation is restricted to the inner few AU of 
the disk, X-rays and FUV photons can operate at tens of AU from the star. 

EUV+X ray and FUV+X-ray models show a similar qualitative b ehavior to 
phot oevaporation by EUV alone, but also several important differences ( Armitag3 
2010l ). Most notably, they predict photoevaporation rates of the order of 10~^ Mq yr" 



which are two orders of magnitude greater than pure EUV photoevaporation 
models. Consequently an inner hole forms early in a disk's accretion history. 
As a result, EUV+X ray and FUV+X-ray models predict a significant popula- 
tion of pre-main sequence stars with relatively massive ( ^ lOMjup) disks with 
large inner holes and no, or very little, accretion. However, observations have 
shown that disks around WTTS tend to have rnuch lower rn asses, ^ 1 — 2M_Jup 
(| Andrews Wilhamd [2OO5I . l2007al . ICieza et al.ll2010l . I2OO8I I. This suggests that, 
in reality, the photoevaporation rates and d isk rn asses are lower th a n tho se pre- 
dicted bv ICorti. Dullemond fc HollenbachI ( 20091 ) and lOwen et al.l ( 20ld ) when 
the inner hole is initially formed and that the quantitative predictions of current 
photoevaporation models should be taken with caution. 
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6.3 Grain growth and dust settling 

Even though sohd particles only represent 1% of the initial mass of the disk, 
understanding their evolution is of utmost interest for disk evolution and planet 
formation studies. Solids not only dominate the opacity of the disk, but also 
provide the raw material from which the terrestrial planets and the cores of the 
giant planets (in the core accretion model) are made. Although viscous accretion 
and the photoevaporation processes discussed above drive the evolution of the 
gas, other processes operate on the solid particles, most importantly, grain growth 
and dust settling. 

Grain growth and dust settling are intimately interconnected processes. Gas 
motions differ slightly from Keplerian motions due to pressure. Small (r ~ 
0.1 /im) grains have a large surface-to- mass ratio and are swept along with the 
gas. As grains collide and stick together, their surface-to. mass ratio decreases 
and their motions decouple from the gas. They therefore suffer a strong drag 
force and settle toward the midplane. This increases the density of dust in the 
interior of the disk, which accelerates grain growth, and results in even larger 
grains settling deeper into the disk. If this process were to continue unimpeded, 
the end result would be a perfectly stratified disk with only small grains in the 
disk surface and large bodies in the midplane. However, because circumstellar 
disks are known to be turbulent, some degree o f vertical stirring and mixing of 



grains is expected (jPullemond Sz Dominikll2005l ). 



6.3.1 Models Grain growth represents the baby steps toward planet for- 
mation. Over 13 orders of magnitude in linear size separate sub-micron particles 
from terrestrial planets, however, and many poorly understood processes oper- 
ate along the way. Idealized dust coagulation models, ignoring fragmentation 
and ra dial drift, predict extremely efficient grain growth. Dullemond &: DominikI 



iiooi) investigated the dust coagulation process in circumstellar disks coupled 



to the settling and turbulent mixing of grains. They included 3 relatively well 
understood dust coagulation mechanisms (Brownian motion, differential settling 
and turbulence) and conclude that these processes are efficient enough to re- 
move all small grains (r < 100 /jui) within 10^ yr. This is clearly inconsistent 
with the wealth of observational evidence showing the presence of micron sized 
grains throughout the duration of the primordial disk phase discussed in the next 
section. The conclusion is that small grains must be replenished and that the per- 
sistence of small gr ains depends on a com plex balance between dust coagulation 



and fragmentation (jPullemond &: Domin ik 2008). 



Recent, more realistic models including fragmentation and radial drift confirm 
the necessity of grain fragmentation to explain the ubiquity of small grains i n disks 
(jBirnstiel. Ormel &: Dullemondll201ll : iBrauer. Dullemond &: Henningjl2008l ) . These 



same models also confirm the severity of the problem known as the "meter-size 
barrier" , a physical scale at which solids are expected to suffer both destructiv e 
collisions and removal through rapid inward migration ( Weidenschillinel 19771 ). 



Even though several possible solutions have be en proposed, including the forma- 



tion of planetesimals in long-lived vortices (e.g.. lHeng &: Kenvoiill2010l ). or via the 
gravi tational instability of millimeter-sized chondrules in a gas deficient subdisk 
(e.g., Youdin fc Shu 20021 ). overcoming this barrier remains on e of the biggest 



challenges for planet formation theories ( Chiang &: Youdin 2010l ) 



6.3.2 Evidence FOR GRAIN GROWTH FROM SUB-MICRON TO MICRONS Pris- 
tine dust in the interstellar medium is composed primarily of amorphous silicates. 
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generally olivine and pyroxene with characteristic features at 9.7 and 18.5 //m 
from Si-0 stretching and 0-Si-O bending modes, respectively. The shape of 
these features is a sensitive diagnostic of the size of small grains, generally in the 
r ~ 0.1 — 5 fim range. Although the emission features of the smallest grains are 
strong and narrow, those of larger particles are weaker and broader. Thermal 
annealing can modify the lattice structure of magnesium-rich olivine and pyrox- 
ene and turn th em into their crystallized forms, known as enstatite and forsterite 
(iHennindlioTnl ). These crystalline silicates have multi-peak features at slightly 
longer wavelengths than their amorphous counterparts. 

Recent surveys with the IRS on Spitzer have allowed the stud y of silicates 
in hundreds of circumstellar disks in nearby star-forming regions ( Fur Ian et al.l 
20061 . l2009l . iKessler-Silacci et al.l[2006l . iMcClure et "al]l2010l . lOliveira et al.ll2010l ). 



Because these observations trace the warm optically thin disk "atmosphere" , and 
the larger particles are expected to settle toward the disk interior, they provide 
information on the smallest population of grains present in the disk. Nevertheless, 
the observed silicate features in most disks are consistent with the presence of 
micron-sized particles and the absence of sub-micron dust grains. This implies 
either that grain growth is more efficient than fragmentation at these scales or 
that sub-micron grains are efficiently removed from the up per layers of the disk 
by stellar winds or radiation pressure (jOlofsson et al.ll2009l l. 

These Spitzer surveys have also shown that the signatures of grain growth and 
crystallization are seen at very early stages of disk evolution, even before the 



envelope has dissipated (iFurlan et al. I l2009l . iMcClure et al]l20inl ). There seems 



to be little connection between the age or evolutionary stage of a primordial disk 
and the dust properties revealed by the silicate features. Even though many 
studies have searched for a correlation between the large scale properties of the 
disk (e.g., mass, accretion rates) and grain characteristics (si zes, degree of crys- 
tallization), no conclusive evidence has yet be en found (e.g., iFurlan et al. 2009 . 
Kessler-Silacci et ahllioO^ . loi iveira et al.ll20ld lV This lack of correlation between 



age and dust properties suggests that the characteristics of the dust population 
depend on a balance between grain growth and destruction and between crystal- 
lization ( via therm al annealing) and amorphization (e.g., via X-ray irradiation, 



Glauser et al]l2009l ^. This balance seems to persist throughout the duration of 



the primordial disk stage, at least in the surface layers of the disk. 
6.3.3 Evidence for grain growth from migrons to millimeters A 
second, independent line of evidence for grain growth is found in the slope, 
ctmm; of the SED at sub-millimeter wavelengths, Fi, oc 'pj^g slope between 

A 0.5 — 1 mm is significantly shallower in protoplanetary disks, Qttitti ~ 2 — 3 
( Andrews Williams 2005 . 2007a . Beckwith &: Sargentlll991 . Mannings fc Emerson 



19941 ) than in the diffuse interstellar medium, amm ~ 4 ( Boulang er et alJ 1199^ 1 



In the Ray leigh- Jeans limit, the flux, Fi, oc B^(l — e 



for optically thick 



emission and — )• Hi^u for optically thin emission. The SED slope is therefore 
bound between amm = 2 and 2 + 13 (see equation [2]) . Unresolved photometric 
observations include some optically thick emission from the inner disk so the con- 
nection between the observed disk average a mm and the grain op acity index, 
requires modeling the surfa ce density profile ( Beckwith et al. 1990) • For a simple 
power law thin disk model, Andrews &: Williamsl ( 20051 ) show that amm — /3 lies 
between 1.1 and 1.8 for disk masses 10~^ — 10~^ Mg^ res pect ively, and is approxi- 
mately independent of a mm • Beckwith &: Sargent ( 199ll ) and Mannings &: Emerson 



Protoplanetary Disks 



31 



(|l994l ) make similar corrections and all reach the same conclusion that the grain 
opacity index is signif icantly lower in dis ks, /3disk ^ 0.5 — 1, than in the interstellar 
medium, /3ism ~ 1-7 (|Li &: Draind 1200 ll ). 

The decrease in the dust opacity index is best explained by the presence of 
substantially larger dust grains in disks relative t o the interstellar medium. For a 
power law d istribution o f graiii sizes, n(a) oc a~^. iD'Alessio. Calvet &: Hartmann 
(|200lh and iNatta et al.l t()04 ) show that (3 ~ 1.7 if the maximum grain size 
flmax < SO/im but then decreases to ^ 1 for Omax ~ 0.5 mm and p = 2.5 — 3.5. 

3)/3isM if flmax is more than 3 
times the observing wavelength, largely independent of grain composition. The 
implication of shallow disk SED slopes at sub-millimeter wavelengths relative to 
the diffuse interstellar medium is that dust grains have grown by at least 3 orders 
of magnitude, from microns to millimeters. 

Some grain growth appears to happen in dense molecular cores before disks 
form. Growth to micron sizes is in fe rred from "co r eshine " , or scattering, of mid- 
infrared light ( Pagani et al. 2010l ). Shirley et al. ( 2000l ) measured the spectral 
slope for 21 cores, some pre-stellar, others in very early stages of low mass star 
formation and find an average spectral slope between 450 and 850 fim, Omm = 2.8, 
which is intermediate between the diffuse interstellar med ium and protoplanetary 
disk s. Multi- wav e length observations of Class cores bv lj0rgensen et al.l ( 20071 ) 
and Kwon et al. ( 20091 ) show opacity indices, /3 1, also shallower than the 
diffuse interstellar medium but not quite as steep as in protoplanetary disks. 

Th ere is also evidence that gra in growth continues throughout the disk life- 
time. Andrews &: Williamsl (|200.5l ) find a significant difference in the SED slope 
between Taurus Class I (median amm = 2.5) and Class II YSO disks (median 
C(mm = 1-8). They also show that the sub- millimeter and infrared SED slopes 
are correlated with a best fit linear relationship, Omm — 

2.09 - 0.40aiR. Confir- 
mation of this single-dish result, which may be biased by cloud contamination in 
the early phases, requires inter ferometry. 

These sub-millimeter data only place a lower bound on the size of the most 
massive grains. Longer wavelength observations constrain the presence of larger 
grains but the flux decreases (as the surface area of the dust decreases for a given 
mass) and th e free-free em ission from an ionized stellar wind often dominates 
for A ^ 1 cm ( Natta et al ][2004). However, the correction for an opaque inner 
region is smalle r . Sur veys of Taurus d isks have b e en car ried out at 7 mm by 
Rodmann et~^ \20()(i ) and at 3 mm by iRicci et al.l I (l2010bl ^ who both find that 



the shallow SED slope generally extends to their observing wavelength and infer 
average opacit y indices, (13) = 1,0.6 respectively. Similar results are found in 
other regions ( Lommen et al.l 2007, Ricci et al. 2011, 2010al ). The implication is 
that centimeter sized particles are commonplace. The longest wave length detec- 



tion o f dust in a protoplanetary disk to date was carried out by IWilner et al. 



( 20051 ). By resolving the nearby, bright TW Hydra disk and showing that its flux 
is constant with time, they convincingly show that the 3.5 cm emission is from 
thermal dust emission and not accretion shocks or a stellar wind. The Draind 
(|2006l ) formulation would therefore indicate the presence of ~ 10 cm sized bodies 
in the disk out at least to the 1 5 AU beam. A similar result has been found for 

WW Cha (iLommen et al.l l2009l'). 

Collisional grain growth ( Blum &: Wurm 20081 ) is expected to have a strong 
radial dependenc e due to the decreasing density a nd rotational velocity with 
increasing radius. Isella. Carpenter fc Sargent ( 2010l ) resolve the RY Tau and DC 
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Tau disks at A = 1.3 and 2.8 mm and model the data with a radially dependent 
dust opacity. However, to within their precision, A/3 = 0.7, they do not see a 
significant gradient. We will soon see detailed images of (3 in disks and learn much 
more about the radial variation of grain growth via multi-wavelength observations 
with ALMA. 

6.3.4 Evidence for dust settling Protoplanetary disks are flared (see 
^4.3.2p and hence intercept and reprocess more stellar radiation than physically 
thin disks. Nevertheless, most T tauri stars exhibit less mid-infrared emission 
than expected for a disk in hydrostatic equilibrium. This can be understood in 
terms of dust settling which r educes the scale height and flaring angle of the disk 



(jPullemond Sz DominikI |2004|) . Mid-infrared slope s can therefore be used as a 



diagnostic of dust settling. D'Alessio et al. ( 20061 ) explored parametric models 



with two populations of grains: micron sized grains with a reduced dust-to-gas 
mass ratio in the surface layers of the disk, and mm-sized grains with an increased 
dust-to-gas mass ratio in the disk midplane. In order to quantify the degree of 
dust settling, they introduced the parameter e, the ratio of the dust-to- gas mass 
ratio in the surface layer of the disk to the standard dust-to-gas mass ratio of the 
interstellar medium (1:100). In this context, e = 1 impli es that no dust settliuj 



has occurred, and e decreases as dust settling increases. iD'Alessio et al.l (120061 ) 
find that the median mid-infrared slopes of CTTS SEDs imply that e < 0.1 and 
p ossibly as low as 10 ~^ 



Fur Ian et al. I (l2006l ) applied these parametric models to reproduce the spectral 



slopes of a sample of over 80 Taurus T Tauri stars observed with the Spitzer 
Infrared Spectrograph. They conclude that most objects are consistent with dust 
depletion factors in the surface layers of the disk of order of 10 to 1000 (i.e.. 



corres ponding to e = 10 — 10 ). A similar result is found by iMcClure et al 



m the Ophiuchus molecular cloud. They also find evidence for significant 



dust settling in young (age ~ 0.3 Myr) objects embedded in the cloud core, 
suggesting that th is process sets in early in the evolution of the disk (see also 
Furlan et al.1l2009l ). 



6.4 Typical evolution and diversity of evolutionary paths 

Although we are far from fully understanding the complex evolution of protoplan- 
etary disks, a coherent picture is starting to emerge from the many models and 
observational constraints discussed above. Although it is clear that not all disks 
follow the same evolutionary path, many observational trends suggest that most 
objects do follow a common sequence of events. In what follows we summarize 
our current understanding of the "typical" evolution of a protoplanetary disk. 

6.4.1 The evolution of a typical disk Early in its evolution, the disk 
loses mass through accretion onto the star and FUV photoevaporation of the 
outer disk ( Gorti. Dullemond &: Hollenbach 20091 ). The FUV photoevaporation 



is likely to truncate the outer edge of the disk, limiting its viscous expansion to a 
finite size of several hundreds of AU in diameter (Figure [6^) . During this "mass 
depletion" stage, which can last several Myr, an object would be classified as a 
CTTS based on the presence of accretion indicators. Accretion may be variable 
on short timescales, but show a declining long-term trend. 

At the same time, grains grow into larger bodies that settle onto the mid-plane 
of the disk where they can grow into rocks, planetesimals and beyond. Accord- 
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ingly, the scale height of the dust decreases and the initiaUy-flared dusty disk 
becomes flatter (Figure [6)3). This steepens the slope of the mid and far-infrared 
SED as a smaller fraction of the stell ar radiation is intercepted by circumstel- 
lar dust ( Dullemond &: Dominik 20051 ). The near- infrared fluxes remain mostly 
unchanged because the inner disk stays optically thick and extends inward to 
the dust sublimation temperature. The most noticeable SED change during this 
stage is seen in the decline of the (sub-)millimeter flu x, which traces the decrease 
in the mass of millimeter and smaller sized particles (j Andrews &: Williamsll2005l . 



2007al . Figure [ZD. 



As disk mass and accretion rate decrease, energetic photons from the stel- 
lar chromosphere are able to penetrate the inner disk and photoevaporation 
becomes important. When the accretion rate drops to the photoevaporation 
rate, the outer disk is no loiiger ab l e to resupply the inner disk with material 
(jAlexander. Clarke &: Pring l3 l2006al : lOwen et al.lEoiol '). At this point, the inner 
disk drains on a viscous timescale ( ^ 10^ yr) and an inner hole of a few AU in 
radius is formed in the disk (Figure [U:) . Once this inner hole has formed, the 
energetic photons impact the inner edge of the disk unimpeded, and the photo- 
evaporation rate increases further, preventing any material from the outer disk 
from flowing into the inner hole. This halts accretion and results in the rapid 
transition between the CTTS and the WTTS stage. 

The formation of the inner hole marks the end of the slow "mass depletion" 
phase and the beginning of the of the rapid "disk dissipation" stage. By the time 
the inner hole is formed, the mass of the outer disk is believed to be ^ 1 — 2 Mjnp , 
as attested bv the low mass of WTTS disks (lAndrews fc WilUamsl I2OO5I . l2007al . 
Cieza et allboiol . [2OO8I ). During this "disk dissipation" stage, the SEDs of the 



WTTS disks present a wid e range of rnorphologies, as expected for disks with 
inner holes of different sizes ( Cieza et al.ll2007l . IPadgett et al.l[200^ . IWahhaj et al. 



Once the remaining gas photoevaporates, the dynamics of the solid par- 



ticles become dominated by radiation effects (Figure Ell). Although the small 
(r ^ 1 iim) grains are quickly blown out by radiation pressure, slightly larger 
ones spiral in due to the Poynting-Robertson effect and eventually evaporate 
when they reach the dust sublimation radius. What is left represents the ini- 
tial conditions of a debris disk: a gas poor disk with large grains, planetesimals 
and/or planets. The fact that the vast majority of WTTS show no evidence for 
a disk (primordial or debris) implies that not every primordial disk evolves into 
a detectable debris disks. Whether some disks remain detectable throughout the 
primordial to debris disk transition or there is always a quiescent period between 
these two stages still remains to be established. 

6.4.2 The age variable Although it is true that there is some discernible 
correlation between the ages of pre-main sequence stars and the evolutionary 
stages of their disks, this correlation is rather weak. In fact, circumstellar disks 
in every stage of evolution, from massive primordial disks to debris disks and 
completely dissipated disks, are seen in st ellar clusters and asso ciations with ages 
ranging from < IMyr to ~ lOMyr (e.g.. iMuzerolle et al.|[201ol l. Diskless WTTS 
in the core of the young Ophiuchus molecular cloud and the relatively old, gas-rich 
TW Hydra are good examples of these extremes. This weak correlation between 
stellar age and disk evolutionary stage can be explained by the combination of 
two factors: the wide range in the duration of the "mass depletion" stage and the 
short time scale of the "disk dissipation" phase. On the one hand, as discussed 
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in ^5.6.H circumprimary disks in medium separation binary systems are likely 
to have truncated disks with small initial masses. Because the "mass depletion" 
stage in such objects is expected to be very short ( ^ 0.3 Myr), their circumstellar 
disks can easily go through all the evolutionary stages described above in less than 
1 Myr. However, initially massive circumstellar disks evolving in isolation can in 
principle remain optically thick at all infrared wavelengths for up to ~ 10 Myr. 
6.4.3 Evidence for alternative evolutionary paths The observed 
properties of most circumstellar disks are consistent with the evolutionary se- 
quence shown in Figure [6j however, some outliers do exist, showing that not 
all disks follow the same sequence of event in an orderly manner. Perhaps the 
most intriguing outliers are accreting objects that have cleared out inner disks 
but retained massive (> lOMjup) outer disks. Examples of such objects, which 
represent a small subgroup of the so-calle d transition disks, inclu de DM Tau, 
GM Aur (|Naiita. Strom Muzeroli^ l2007l ) . and RX J1633. 9-2442 (ICieza et al. 



2010l ). Their large accretion rates, inner radii, and disk masses make their holes 



incompatible with the evolutionary sequence discussed above. We must then 
consider an additional agent, most likely dynamical clearing by a (sub)stellar 
or planetary-mass companion. The nature of these objects is discussed in some 
detail in the following section. 



6.5 Section summary 

• Protoplanetary disks evolve through a variety of processes, including viscous 
transport, photoevaporation from the central star, grain growth and dust 
settling, and dynamical interaction with (sub)stellar and planetary-mass 
companions. 

• Photoevaporative flows from disk surfaces have been observed but the mod- 
els disagree on the relative importance of FUV, EUV, and X-ray photoe- 
vaporation. 

• There is strong evidence for grain growth to millimeter (and, in some cases, 
centimeter) sizes but the presence and distribution of larger bodies remain 
unconstrained. 

• Most protostellar disks go through a slow "mass depletion" phase followed 
by a rapid "disk dissipation" stage. As the accretion rate steadily drops 
below the photoevaporative rate, the disk is rapidly eroded from the inside- 
out. A wide range in the duration of the two phases, together with the 
intrinsic dispersion of disk masses and radii, weakens the correlation be- 
tween stellar age and disk evolutionary stage. 



7 TRANSITION DISKS 



Transition disks were first identified by IRAS as objects with little or no excess 



emission at A < lOmn and a significant excess at A > 10 /im (jStrom et al.l 11989 



Wolk fc Waited Il996l ). The lack of near-infrared excess was interpreted as a 



diagnostic of inner disk clearing possibly connected to the early stages of planet 
formation. Because of this possible connection, transition disks have received 
special attention in circumstellar disks studies even though they represent a small 
percentage of the disk population in nearby star-forming regions. In this section 
we discuss the diversity of SED morphologies presented by transition disks, their 
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incidence, and their connection to planet formation and other disk evolution 
processes. 



7.1 SED diversity and interpretation 

Transition disks present a variety of infrared SED morphologies and their diversity 
is not properly captured by the traditional classification scheme of young stellar 
objects (i.e., the Class I through III system, s ee !^[2]). In order to better describe 



the shape of transition disk SEDs. lCieza et al.l (120071 ) introduced a two-parameter 



scheme based on the longest wavelength at which the observed flux is dominated 
by the stellar photosphere, Aturn-off, and the slope of the infrared excess, Oexcess, 
computed from Atum-off to 24 ^m (Figure [8^.) 

Because of their diversity, the precise definition of what constitutes a transition 
object found in the disk-evolution literature is far from homogeneous. Transition 
disks have been defined as objects with no detectable near -infrared excess, steep 
slopes in the mid-infrared, and large far-inf rared excesses (jMuzerolle et al.ll2010l . 
Sicilia-Aguilar. Henning fc HartmannI 2010l ) . This definition is the most restric- 
tive one and effectively corresponds to objects with Aturn-off ~ 4.5 — 8.0 /um and 
CKexcess > 0. The above definition has been relaxed by some researchers (Brown 
et al. 2007; Merin et al. 2010) to include objects with small, but still detectable, 
near-infrared excesses. Transition disks have also been more broadly defined in 
terms of a significant decre ment relative to the median SED of CTTS at an y 
or all infrared wavelengths ( Cieza et al. 2010 : Naiita. Strom &: Muzerolle 20071 ). 
Several nouns h ave recently emerg ed in the literature to distinct sub-classes of 
transition disks ( Evans et al. I l2009al ). Objects with no detectable near-infrared ex- 
cess and Oexcess > have been called "classical" transition disks ((Muzerolle et al. 



2010l ). whereas other obje cts with a sharp r i se in their raid-iri frared SEDs have 
been termed "cold disks" ( Brown et al. 2007 . Merm et al. 20 id ) regardless of the 
presence of near-infrared excess. Disks with a significant flux decrement at all in- 
frared wavelengths relative to the SED of an optically thick disk extending to the 
dust sublimatior i temperature typ ically have Oexccss < and been r eferred to as 
disks (ILada et al.l 120061). "homologously dep leted" disks (ICurrie et al. 



anem ic aisKs ll.Lada et ai.Mzuupii. nomoiogousiy aep ietea aisKs in^urrie et al. 
2009l V or "weak excess" disks (jMuzerolle et al.l l2010l ) Finally, disks with evi- 



dence for an optically thin (to starlight) gap separating optically thick inner and 



outer disk components have been called "pre-transition" disks (lEspaillat et al 



20071 . l20inl ;i because they are believed to be precursors of objects with sharp, but 
empty inner holes (i.e., the "classical" transition disks mentioned above). 

The inherent diversity of transition disk SEDs is encapsulated in Figure Eb 
which shows how the range of slopes, aexcessj increases as the point at which excess 
emission is detected, Aturn-off, increases. Disks without inner holes, Aturn-off ~ 2 /im, 
have similar mid-infrared SED slopes but this is not the case for larger Aturn-off 
suggesting a range of processes by which the inner hole is formed and develops. 

Examples of the different types of transition disk SEDs are shown in Figure [H 
The SEDs that decrease in the mid-infrared, aexcess ~ 0, have generally been in- 
terpreted as evidence for grain growth and dust settling toward the midplan e 
(|Cieza et al.ll2010l . iLada et al.l[2006l . ISicilia-Aguilar. Henning &: Hartmannlboid ). 
Both processes tend to reduce the mid-infrared fluxes. Grain growth removes 
small grains and thus reduces the opacity of the inner disk. Dust settling re- 
sults in flatter disks that intercept and reprocess a smaller fraction of the stellar 
radiation. This group of objects includes anemic, homologously depleted disks. 
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and weak-excess disks. Objects with rising mid-infrared SEDs, Oexcess > 0, a 
group that includes classical transition disks, cold disks, and pre-transition disks, 
are more consistent with a steep radial dependence in the dust opacity resulting 
in a sharp boundary to the inner opacity hole. These may be due to dynami- 
cal clearing by a co mpanion, possibly planetary ( Artvinowicz Sz Lubowl 1994 ) , or 
photoevap oration (lAlexander. Clarke &: Pring Debris disks also have 

sharp inner holes but they can be distinguished by their low luminosities and lack 
of gas. 



7.2 Incidence 

Establishing the incidence of transition disks is not straightforward given the dif- 
ferent definitions used by different studies. Also, samples are usually not complete 
and thus suffer from selection effects. Furthermore, several background objects 
can mimic the SEDs of some types of transition disks and skew the statistics. In 
particular, Asymptotic Giant Branch stars an d classical Be star s can be easily 



confu sed with "weak excess" transition disks (jCieza et al.l |2010j, lOliveira et al. 



20101 ) ■ and the SEDs of edge-on protoplanetary disks can look like those of cold 



disks (|Merm et al.ll2010l ). 



With these caveats, the fraction of disks in nearby star-fo rming regions that 



are seen in a transition stage is thought to be ^ 10 — 20% ( Currie et al. 2009 



Dahm fc CarDenterll2009l . iFang et al.l[2009l . iHernandez et al.ll2007l . lKim et al.ll2009l . 
Lada et al. 20061 ). Also, it is clear that "weak-excess" (or "anemic" and "homol- 
ogously depleted" ) transition disks outnumber, by factors o f ~ 2 — 3, objects with 
sharp inner holes (jCieza et al.ll20ld . iMuzerolle et al.ll2oTol '). The fraction of cold 
disks in nea rby star-forining regions has been rece ntly established to lie between 
5 and 10% (|Merm et al.ll2010l . lOliveira et al.ll2010l ). 



Weak-excess disks seem to be more common in older clusters (Muzerolle et al 



2O10l V This is expected because disks should become flatter and lose mass with 
time (see §6.4p . Weak-excess disks appear to be much more abundant around 
young (age ~ 1 Myr) M-s tars than around solar-type stars of the same age 



(jSicilia-Aguilar et al.l 



20081 ). However, this may be a luminosity effect rather 
than an evolutionary one. Because M-type stars are much fainter and cooler 
than solar-type stars, they may present weak mid-i nfrared excess emission even if 
the d isk extends in to the dust sublimation radius ( Ercolano. Clarke fc Robitailld 
20091 ). 

The relatively small number of objects seen in a transition stage suggests that 
the evolutionary path through a transitional disk is either uncommon or rapid. 
However, observations show that an infrared excess at a given wavelength is al- 
ways accompanied by a larger excess at longer wavelengths, out to ~ 100 ^m. This 
implies that, unless some disks manage to lose the near-, mid-, and far-infrared 
excess at exactly the same time, the near-infrared excess always dissipates before 
the mid-infrared and far-infrared excess do. No known process can remove the 
circumstellar dust at all radii simultaneously, and even if grain growth or dynam- 
ical clearing do not produce an inner hole, photoevap oration by the central star 
will do so once the accretion rate through the disk falls below the photoevapora- 
tion rate. Therefore, it is reasonable to conclude that transition disks represent 
a common but relatively short phase in the evolution of a circumstellar disk. Of 
course, how long the transition stage lasts for a given disk strongly depends on 
the nature of each object. On the one hand, the transition phase should be very 
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short, < 0.5 Myr, if the inner hole is formed through photoevaporation. However, 
an object could in principle show a transition disk SED for a longer period of 
time if the inner opacity hole is due to grain growth, or giant planet formation. 



7.3 Physical properties 

Although it is clear that many different processes can result in transition disk 
SEDs, the relative importance of these processes is not well understood. Dis- 
tinguishing among processes requires additional observational constraints on the 
physical properties of these objects, such as disk mass, accretion rates, fractional 
disk luminosities (Ldisk/-^^star), and multiplicity information. These observational 
constraints have only recently been obtained for statistically significant samples 

of transition disks. 

Naiita. Strom &: Muzeroll^ ( 200?! ) investigated the SEDs, disk masses, and ac- 



cretion rates of over 60 Taurus pre-main sequence stars and identified 12 transi- 
tion disks. Based on their SEDs, 7 of these 12 could be further classified as weak 
excess disks, 3 as classical transition disks, and 2 as pre-transition disks. They 
found that the transition disks in their sample have larger average disk masses 
and lower average accretion rates than non-transition disks around single stars. 
They concluded that, with the exception of CoKu Tau 4, the SEDs were all more 
consistent with the giant planet formation scenario than with photoevaporation 
models (which at the time predicted negligible accretion rates and lower disks 
masses) and grain growth (which would favor higher accretion rates for a given 
disk mass). However, because their sa mple was drawn fro m the Spitzer spectro- 
scopic survey of Taurus presented by Fur Ian et al. ( 20061 ). who in turn selected 



their targets based on mid-infrared colors from IRAS, their sample was biased 
towards the brightest objects in the mid- infrared and are unlikely to represent 
the overall population of transition di s ks. 

As a counter example, Cieza et al. (j2008l ) studied over 20 WTTS with tran- 



sition disk SEDs, found that all had very small disk masses ( ^ 2Mjup), and 

concluded they were all con s istent wit h the EUV photoevaporation models of 

Alexander. Clarke fc Pringle (j2006al lbl) . ISicilia-Aguilar. Henning fc Hartmann (l2010l ^ 



studied accretion rates in a sample of 95 members of the ~ 4 Myr Tr 37 cluster. 
They found that half of the 20 classical transition disks in their sample had evi- 
dence for accretion (i.e. their holes contain gas but no dust) and half were non- 
accreting (i.e., their holes are really empty). Furthermore, the accretors had rates 
that were indistinguishable from those of regular C TTS in the cluster, a result 



that i s at odds with the Taurus sample studied by iNajita. Strom &: Muzerolle 
()2007l ). These discrepancies show that transition disks are a highly heteroge- 
neous group of objects and that the mean properties of a given sample are highly 
dependent on the details of the sample selection criteria and, thus, should be 
in terpreted wi t h cau tion 



Cieza et al. I (|2010l l studied a sample of 26 transition disks in Ophiuchus. With 



the exception of edge-on disks, highly embedded objects, and the lowest mass 
transition objects, all of which are too faint in the optical to pass one of their 
selection criteria, the sample is likely to be representative of the entire transi- 
tion disk population in the Ophiuchus molecular cloud. They find that 9 of the 
26 targets have low disk masses (< 2.5Mjup) and negligible accretion, and are 
thus consistent with photoevaporation. 4 of these 9 non-accreting objects have 
fractional disk luminosities < 10"'^, however, and could already be in a debris 
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disk stage. The other 17 objects in the sample are accreting. 13 of these are 
consistent with dust setthng and grain growth (oexcess ~ 0), whereas the other 4 
have rising mid-infrared SEDs (aexccss > 0) indicative of sharp inner holes, and 
are candidates for harboring embedded giant planets. A decision tree showing 
the different possibilities is shown in Figure [TOl 

Interestingly, non-accreting objects with relatively massive disks ( ^ 2.5 Mjup) 
are absent in the Ophiuchus transition disk sample. Such objects are rare bu t 
do exist. The "cold disk" around T Cha is an example ( Brown et al.l 20071 ). 



and its properties could be explained by a companion massive enough to disrupt 
accretion onto the star. A few accreting binary sy stems are kno wn, however, for 



3ys 

exam ple, DQ Tau (jCarr. Mathieu k Najital l200ll ) and CS Cha (jEspaillat et al. 



20071 ). which implies that a massive companion does not always stop accretion. 



The eccentricity of the orbit, and the viscosity and scale height of the disk also 



deter mine whether accretion onto the star can continue (jArtvmowicz Sz Lubow 
19961 ). 



Binarity is also known to result in transition disk SEDs. Identifying close 
companions is particularly difficult, however, especially if they are within the 
inner ~ 5 — 10 AU of the disk, which is beyond the near-infrared diffraction 
limit of 8 — 10 m class telescopes at the distance of the nearest star-forming 
regions (~ 125 — 140 pc). Using t he aperture masking technique in a single Keck 
telescope, Ireland Kraud ( 20081 ) showed that CoKu Tau 4 is a near-equal mass 



binary system, which explains the inner hole that had been inferred from its 
SED; nevertheless, it is clear that not all sharp inner holes are due to binarity. 
DM Tau, GM Aur, LkCa 15, UX Tau, an d RY Tau have all been observed with 
the Keck interferometer ( Pott et al. 2010l ). For these objects, all of which show 



evidence for accretion unlike CoKu Tau 4, stellar companions with flux ratios 
< 20 can be ruled out down to sub-AU separations. 



7.4 Resolved observations 

Although SED modeling provides strong evidence for the presence of inner holes 
in some transition disks, this evidence is indirect and model-dependent. Long 
baseline interferometry at (sub-)millimeter wavelengths, however, can directly 
image the inner holes of some transition disks in nearby star- forming regions. 

Inner holes have finally been resolved in a number of transition disks in recent 
years. Figure [TT] displays a montage of images at the same linear scale in com- 
parison to the orbits of the giant planets in the Solar System. The holes can be 
clearly seen in the face-on disks and are apparent through a central dip in the 
emission for edge-on disks. 

Because the emission is optically thin in these (sub-)millimeter data, the im- 
ages directly show the mass surface density profile of the disk and can be used to 
measure the size and sharpness of inner holes. In most cases, the (sub-)millimeter 
interferometry confirms the expected sizes of the inner holes (ranging from ~ 4 AU 
to ~ 50 AU in radius) and their sharpness (i.e., a large increase in the dust density 



over a sma ll range in radii) predi cted from the SED. For example, iBrown et al. 



(|2009l ) and [Andrews et al. I (|2010bl ) show that the flux decrement toward the cen- 



ters of their imaged transition disks are consistent with step functions with dust 
depletion factors ranging from 10^ to 10^. 

Several of these transition disks are gas rich. In two cases, GM Aur and J1604- 
2130, where the millimeter CO emission is resolved to sufficient detail either 
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spectroscopically or spatially, the cen tral cavity is found to b e not only dust but 
also g as deficient ( Dutrev et al. 20081 . respectively). However, Pontoppidan et al. 
( 20081 ) found mid-infrared CO emission from well inside the dust gaps in SR 21, 
HD 135344, and TW Hya. 

The details in millimeter interferometric images may contain dynamical hints 
of giant planets. All the images in Figure [TT] appear to be asymmetric but these 
are no more than suggestions given the modest sign al-to-noise ratio and poor 
sampling of the Fourier plane. iHughes et al. note a warp between the 

dust and gas structures in GM Aur. High fidelity images with ALMA will soon 
allow us to study the structure of circumstellar disks in much greater detail, and 
the incidence of radial gaps, density waves, warps, and other possible azimuthal 
asymmetries should become evident. It may well be possible to find earlier stages 
in the disk clearing process throug h a decline in the m illimeter flux before the 
infrared opacity drops below unity ( Andrews et al.l 20081 ). 

In addition to the (sub-)millimeter images discussed above, in at least one 
case, LkCa 15, the inner hole has been imaged in the near -infrared from starlight 
reflected from the edge of a disk wall ~ 50 AU in radius ( Thalmann et al. 2010l ) 
once again confirming the structure inferred from modeling the SED of the object. 



7.5 Section summary 

• Transition disks can be broadly defined as disks with a significant flux 
decrement relative to the median SED of CTTS at any or all infrared wave- 
lengths. They constitute at most 20% of the disk population. 

• There is a wide range of transition disk SEDs, indicative of the varied phys- 
ical processes, photoevaporation, grain growth, and dynamical interactions 
with companions or planets, that produce them. 

• Grain growth and dust settling produce SEDs with falling mid-infrared 
emission (oexcess ~ 

0). SEDs that rise in the mid-infrared (ctexcess > 0) &re 
more consistent with a sharp boundary to the inner hole due to photoevap- 
oration or dynamical interactions. 

• Giant planet formation best explains the combination of accretion and 
steeply rising mid- infrared emission in moderately massive disks, M ^ 3 Mjup. 



8 SUMMARY POINTS 



Observations of protoplanetary disks are challenging due to their small size, low 
masses, and cool temperatures. However a number of basic facts have been clearly 
established. Mid-infrared observations of optically thick emission provide the 
most sensitive measures of the presence of a disk through which we infer their 
occurrence and lifetimes. Millimeter wavelength observations of optically thin 
emission provide the best measures of disk masses and, through interferometry, 
resolved images of their structure. The interpretation of the data is complicated 
by the uncertainties in our knowledge of grain growth and settling, gas dispersal, 
and the feedback between composition and structure. Nevertheless, it is clear 
that there are a variety of evolutionary pathways, and many different physical 
processes competing with each other, including planet formation. Even the re- 
stricted subject area of this review on observations of the outer disk is so large 
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that we have given a summary to each section. Here, we distill these summaries 
yet further into six main points. 

• Circumstellar disks form almost immediately after a molecular core col- 
lapses. They appear to be highly unstable at early times and accrete in 
bursts onto the central protostar. 

• By the time a YSO becomes optically visible, the mass of the now pro- 
toplanetary disk averages 1% that of the central star. Its surface density 
increases approximately inversely with radius and it has a soft, exponen- 
tially tapered edge between a few tens to hundreds of AU. The disks are 
now generally stable and rotate with a Keplerian velocity profile. About 
15% of disks around solar mass stars have a MMSN (lOMjup) of material 
within 50 AU radius. 

• Disks around solar and lower mass stars have a median lifetime between 2 
and 3 Myr but with a large dispersion from less than 1 Myr to a maximum 
of 10 Myr. Lifetimes are shorter around higher mass stars and binaries 
with semi-major axes between 5 and 100 AU. The large scale environment 
is relatively unimportant: disruption by stellar flybys is very rare and pho- 
toevaporation by massive stars generally only erodes the outer disk, beyond 
about 50 AU. 

• Protoplanetary disks evolve through a variety of processes, including viscous 

transport, photoevaporation by the central star, grain growth and dust 
settling, and dynamical interaction with (sub)stellar and planetary-mass 
companions. Most disks evolve via a slow decrease in the mass of gas and 
small particles followed by rapid disk clearing at all radii. 

• The growth of dust grains from sub-micron sizes in the ISM to millimeter 

sizes in disks occurs early and continues with time. The presence of snow- 
balls or pebbles several centimeters in size is hard to measure but has been 
inferred in a couple of cases. 

• About 10-20% of disks show mid-infrared dips indicative of inner holes. 

They present a wide range of SED morphologies and physical properties 
(disk masses, accretion rates, and inner hole sizes). These transition disks 
are excellent laboratories to study the disk clearing phase and the formation 
of planets. 

9 FUTURE ISSUES 

Protoplanetary disks emit predominantly at long wavelengths where the atmo- 
spheric background is high, and where high resolution requires interferometry. 
By necessity, many current observations are of the brightest disks but these may 
not be truly representative. Systematic studies of the median disk population are 
about to becomre routine. Ongoing observations with Herschel at far-infrared 
wavelengths around the peak of their SED are just beginning to produce new 
insights into disk structure, chemistry, and evolution. Within the next few years, 
ALMA will dramatically expand upon the pathfinder work of present-day interfer- 
ometers at (sub-)millimeter wavelengths. Further down the line lies the prospect 
of sensitive, high resolution observations in the mid-infrared with thirty-meter 
class telescopes, as well as at centimeter and longer wavelengths with the Square 
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Kilometer Array. There is an enormous range of issues to investigate with these 
and other instruments. Here, we hst those that we find particularly compelling. 

• Disk formation: With the high sensitivity and imaging fidelity of ALMA, 
it will be possible to map faint isotopic lines of dense gas tracers and search 
for small rotationally supported structures in the centers of molecular cores. 
Instabilities may be revealed through spiral waves and other asymmetries 
and can provide an independent dynamical disk mass estimate. Optically 
thick lines seen in absorption against the disk continuum will show the infall 
of material from the core. Multi-wavelength continuum imaging will track 
the increase in the grain size distribution through the process. 

• Peering into the terrestrial planet zone: Very high resolution observa- 
tions, ^0.1, at wavelengths beyond a millimeter can image optically thin 
dust emission and resolve structures in the terrestrial planet-forming zone, 
R ^ 5 AU, of the closest disks. The longer the wavelength of the observa- 
tion, the larger the size of the dust grains that can be detected. The first 
systematic studies at centimeter wavelengths arc just beginning with the 
extended VLA. The Square Kilometer Array will measure the distribution 
of rocks and snowballs up to meter sizes. Together with ALMA, resolved 
images from sub-millimeter to centimeter wavelengths will show the radial 
variation of grain growth. 

• Disk chemistry: The ability to survey many disks in many lines with ALMA 
will revolutionize the young field of disk chemistry. As different species and 
transitions are excited in different regions of a disk, such observations will 
enable a far more complete picture of the gas disk structure to be developed. 
Observations of H2^0 will reveal the water context of disks and constrain 
the location of the snowline in a statistically meaningful sample. The detec- 
tion of other molecular isotopologTics will allow isotopic abundances (and 
radial gradients) to be measured and directly compared to cosmochemical 
studies of meteorites. 

• What is the overall evolution of the gas-to-dust mass ratio in circumstel- 
lar disks? Although the evolution of the dust (r ^ 1mm) content in disks 
can be traced reasonably well by current (sub)millimeter observations, our 
current understanding of the evolution of the gas is very limited. Obser- 
vations of gas tracers with Herschel (e.g., [O I] at 63.2 /um), large infrared 
telescopes (e.g., H2 at 12.4 and 17.0 /um), and ALMA (e.g., rotational lines 
of CO and its isotopologues) will reveal the evolution of the gas and the 
gas-to-dust ratio, which is critical to understanding the formation of both 
terrestrial and giant planets. 

• How and when do giant planets form? Although much progress has been 

made in understanding the structure and evolution of circumstellar disks, 
this fundamental question still remains unanswered. Detailed studies of 
embedded (Class I) disks with ALMA will help to establish whether massive 
young disks can be conducive to the formation of giant planets through 
gravitational instability. Similarly, planet searches in Class II YSOs, and 
transition disks in particular, with thirty-meter class infrared telescopes will 
identify young giant planets at the last stages of the core accretion process. 

• Toward comprehensive disk evolution models: To date most disk evolution 
models have focused on one or two physical processes at a time (e.g., viscous 
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accretion and photoevaporation) while ignoring other equahty important 
ones (e.g., grain growth and dust settling and dynamical interactions with 
(sub)stellar companions and/or young planets). In reality, however, it is 
clear that all these processes are likely to operate simultaneously and affect 
one another, and that any realistic disk evolution model should include all 
known disk evolution mechanisms. 

• Placing our Solar System in context: The small number of protoplanetary 
disks that have been studied in detail each have their own idiosyncracies. 
Whereas there are rough matches to the mass, size, and surface density 
profile of the MMSN, it is not clear how common these conditions were. Im- 
provements in technology at all wavelengths will allow more refined studies 
of the closest disks and large surveys of more distant star-forming regions. 
In tandem with the ever increasing knowledge-base on the number, mass, 
and density distribution of cxoplancts, wc will gather the detailed statis- 
tics necessary to understand what is typical and what is atypical about the 
protosolar nebula. 
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Table 1: Classification of Young Stellar Objects 



Class 


SED slope 


Pliysical properties 


Observational characteristics 







Menv > Mstar > Mdisk 


no optical or near-infrared emission 


I 


aiR > 0.3 


Mstar > Menv ~ Mdisk 


generally optically obscured 


FS 


-0.3 < aiR < 0.3 




intermediate between Class I and II 


II 


-1.6 < aiR < -0.3 


Mdisk/Mstar ~ l%,Afc„v ~0 


accreting disk; strong Ha and UV 


III 


am < -1-6 


Mdisk/Mstar <. 1%, Menv ~ 


passive disk; no or very weak accretion 




Figure 1: The distribution of protoplanetary (Class II) disk masses in Tau- 
rus, Ophiuchus, and Orion. The dust masses are derived from milhmeter fluxes 
and extrapolated to a total mass assuming a maximum grain size of 1 mm, 
characteristic temperature of 20 K, and an interstellar ga s-to-dust ratio of 100 
([Andrews Williamsl[2005l . l2007al . iMann fc Williamsll2010l ) . The filled bars show 
the range where the millimeter measurements are complete for each region. The 
hashed bars in the Orion histogram show the disks with projected distances 0.3 pc 
beyond the 06 star, 0^ OriC, in the Trapezium Cluster. 
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Figure 2: Radial surface density (gas+dust) profiles for Class II YSO disks in 
Ophiuchus based on fitting an exponentially tapere d po wer law profile t o 880 fim 
visibilities and infrared SEDs ( Andrews et al. 20091 ) and Andrews et al. I (12010^ ). 
The dark gray rectangular regions mark the MMSN surface densities for Saturn, 
Uranus, and Neptune. 



Protoplanetary Disks 



57 



T — ■ — ' — ■ — I — ■ — ' — ■ — I — ■ — ' — ■ — I — ■ — ' — ■ — r 




J 1 1 1 I 1 1 1 I 1 1 1 I 1 1 1 L 



4 2 -2 -4 4 2 -2 -4 

Aa (■■) Aa (") 

Figure 3: CO(3-2) emission from the disks around TW Hydra (left) and 
HP 163296 (right) observed with the SMA at a spectral resolution of 44ms~^ 
( Hughes et al.ll2010l ). The contours show the zeroth moment (velocity-integrated 
intensity), whereas the colors show the first moment (intensity- weighted veloc- 
ity). The synthesized beams are shown in the lower left corner of each panel with 
a size of 1. 7 x 1. 3 at position angles of 19° and 46° degrees for TW Hydra and 
HD 163296, respectively. The contours start at 3a and increase by intervals of 
2a, where the rms noise a = 0.6 Jybeam"^. 
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Figure 4: Sub-millimeter spectroscopy of molecular rotational lines in the chem- 
ically rich nearby TW Hydra disk. These observations, made with the SMA are 
at a range of resolutions, shown in the lower left corner of each panel. The nearly 
face-on disk generally shows centrally peaked emission except for the DCO"*" line 
which peaks in a ring where the temperature is colder than the inner disk and 
CO freezes out of the gas phase onto grain surfaces. (Figure courtesy of Charlie 
Qi.) 
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Figure 5: The variation of protoplanetary disk mass with the mass of the central 
star. Upper hmits are only shown at the extremes of the stellar mass range where 
no disks have been detected. The dashed diagonal lines delineates where the mass 
ratio is 1%, and is close to the median of the detections. Almost all the disks 
around stars with masses M* = 0.04 — 10 Mq lie within the grey shaded area, 
±1 dex, about the median. The exception are O stars where no disks are detected 
at (sub-)millimeter wavelengths, indicating either very short disk lifetimes or a 
different star formation scenario. 
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Evaporation flow 




Evaporation flow 



c) Photoevaporating disk d) Debris disk 




Figure 6: The evolution of a typical disk. The gas distribution is shown in blue 
and the dust in brown, (a) Early in its evolution, the disk loses mass through 
accretion onto the star and FUV photoevaporation of the outer disk, (b) At the 
same time, grains grow into larger bodies that settle to the mid-plane of the disk, 
(c) As the disk mass and accretion rate decrease, EUV-induced photoevaporation 
becomes important, the outer disk is no longer able to resupply the inner disk with 
material, and the inner disk drains on a viscous timescale (~ 10^ yr). An inner 
hole is formed, accretion onto the star ceases, and the disk quickly dissipates from 
the inside out. (d) Once the remaining gas photoevaporates, the small grains are 
removed by radiation pressure and Poynting-Robertson drag. Only large grains, 
planetesimals, and/or planets are left This debris disk is very low mass and is 
not always detectable. 
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Figure 7: Evolution of the disk masses across the empirical sequence defined 
by the slope of the infrared SED. The circles show the median disk mass, as 
measured at sub-millimeter wavelengths for a sample of 300 YSOs in Taurus 
and Ophiuchus. Error bars show the distribution quartiles. Transition disks are 
defined here as those YSOs that lack infrared excesses for wavelengths less than 
25 /im, but have detectable sub-millimeter emission. With this definition, no 
Class III YSO was detected and the string ent limit to their median ma ss comes 
from stacking the non-detections together ( Andrews &: Williams 2007al ) . 
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Figure 8: (a) A schematic of the Oexcess and Aturn-off classification for a tran- 
sition disk. Although the aiR value is typical of a CTTS with a full disk, the 
long Aturn-off wavelength and positive aexcess indicate the presence of an inner 
hole, (b) Distribution o f «pyrsRs with respec t to Aturn-off for a range of different 
disk evolutionary types ( Harvey et al. 200?! ) . The diagram demonstrates the in- 
creasing range of possible Oexcess values at longer Aturn-off wavelengths due to the 
diversity of transition disk SED morphologies. 
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Figure 9: The diversity of transition disks SEDs. The points represent pho- 
tometry from optical to mid-infrared wavelengths, whereas the dark solid lines 
are Spitzer infrared spectra. The stellar photosphere is shown as a blue curve, 
and the dark hashed region shows the range of SEDs for typical accreting T 
Tauri stars, (a) A weak-excess, anemic, or homologously depleted disk has a 
significant flux decrement at all mid-infrared wavelengths relative to the T Tauri 
SED. (b) A cold disk or "classical" transition disk displays excess emission above 
the photosphere only at mid-infrared wavelengths and beyond, (c) A cold disk 
with little near-infrared emission and a strong 10 /im silicate feature, (d) A 
cold disk with near-infrared excess emission. This can also be considered a pre- 
transition disk because its SED can be modeled with an optically thin gap sep- 
arating optically thick inner and outer disk components. Figure adapted from 
Najita. Strom fc Muzeroli^ (|2nn7l ). 
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Figure 10: Decision tree to determine the dominant physical process for transi- 
tion disks, defined as having a significant flux decrement relative to the median 
SED of CTTS at any or all infrared wavelengths. This ap pears to be an exh aus- 
tive set of possibilities for transition disks in Ophiuchus ( Cieza et al. 20 id ) but 
there may be rare additional possibilities in other regions. 
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Figure 11: Montage of (sub-)millimeter images of protoplanetary disks with 
resolved inner holes. Each panel is 400 AU on a side, and the images are shown 
stretched over 99% of the range from minimum to maximum. The gray ellipse 
at the lower left of each panel shows the synthesized beam size, and ranges from 
0. 15 to 0. 8. The orbits of the giant planets and Pluto in the Solar System 
are shown in the lower right panel for scale. The data are from a variety of 
published studies with the SMA, Plateau d e Bure, and CARMA interferometers 



(lAndrews et~aD l2010bl. iBrown et al.l l2009l. iGuilloteau et all l2008l . iHughes et al 



2009a, Isella. Carpenter fc Sargent 2009, Pietu et al. 20061 ) as well as not yet 



published results. 



